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Introduction 

Getting  Started/Operation 


In  the  late  1980s  S.  Jenson  from  the  Eros  Data  Center  developed 
a  series  of  FORTRAN  programs  to  extract  surface  hydrologic 
information  from  digital  terrain  data  (Jenson  and  Domingue, 
1988).  The  Hydrologic  Terrain  Analysis  Software  (HTAS)  adds 
to  the  existing  GIS  terrain  analysis  capabilities  a  series  of  tools 
used  to  aid  sound  design  of  resource  and  environmental  plan- 
ning and  management  projects  which  utilize  information  about 
surface  hydrologic  processes.  Among  its  many  capabilities, 
HTAS  permits  discovery  of  surface  depressions,  flow  accumula- 
tion, flow  directions,  and  drainage  pattern  and  watershed 
delineation,  with  a  choice  of  manual  or  automatic  selection  of 
potential  watershed  outlets.  Potential  applications  include  a 
variety  of  modeling  goals,  such  as,  runoff  and  sediment  yield 
modeling,  stream-flow  modeling,  dam  and  reservoir  design, 
erosion  control  structure  design  (such  as  spreaders),  riparian 
area  studies  including  input  of  sediment  and  pollutants  into 
streams,  and  aquifer  recharge  area  discovery.  Some  of  these 
applications  have  been  or  are  being  developed,  and  others  are 
yet  to  be  identified. 

This  manual  provides  an  introduction  to  the  concept  and  the 
capabilities  of  the  HTAS,  as  well  as  step  by  step  instructions  on 
its  use.  Although  at  this  time  the  BLM's  Map  Overlay  and 
Statistical  System  (MOSS)  and  Map  Analysis  Processing  System 
(MAPS)  software  do  not  include  the  HTAS  capabilities  as  com- 
mands, and  the  GIS  interface  remains  cumbersome,  we  believe 
that  the  benefits  of  using  the  software  offset  the  difficulties. 

Program  Installation  and  Execution 

The  HTAS  consists  of  14  FORTRAN  programs  which  run  under 
the  Primos,  UNIX,  or  VMS  operating  systems.  Each  program 
compiles  and  runs  as  an  independent  unit,  consequently  no 
subroutine  libraries  are  required.  The  attached  diskette  contains 
the  executable  version  of  the  programs  compiled  under  Primos. 
Create  a  directory  to  store  these  programs  and  copy  them  off  of 
the  diskette.  To  execute  the  programs  use  the  Primos  RESUME 
(R)  command  using  the  format: 

R  (path)  (f  i  lenatne) 

I  If  starting  at  the  directory  where  you  stored  the  programs  e.g., 
^  IS>HVDRO,  type: 

R  CONCRUE 
If  starting  at  any  other  location  type  the  entire  pathname: 

R   IS>HVDRO>CONCRUE 
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Prompting 

The  programs  prompt  for  various  parameters,  including  the  size 
of  the  map  being  processed  in  lines  (rows)  and  samples  (col- 
umns), file  names  and  system  name  before  processing  begins. 
Enter  the  file  names  and  system  name  (PRIME)  enclosed  in 
single  quotes,  e.g.,  'PRIME'. 

R  DELTR 

ENTER  NL,NS,DIR  FILE, COUNT  FILE, DELTA 
FILE, SVSTEn<UNIX,UMS, PRIME) 

1334,  1278,  'DIR.HflP',  'LEU.nflP',  'DELTA. MflP'/PRiriE' 

File  cliaractepistics  and  fiia  transfer 

The  programs  operate  on  binary  1*2  and  1*4  disc  files.  They 
require  direct  access,  fixed  length  record  file  type  to  access 
information  at  any  place  in  the  file  at  will.  The  program  CON- 
VERT permits  transfer  of  IDIMS  files  into  direct  access,  fixed 
length  format.  The  MOSS/MAPS  interface  part  contains  addi- 
tional information  on  this  process. 

All  the  raster  data  sets  exist  in  the  1*2  format  (cell  values  are 
represented  as  an  integer  value  stored  in  a  two  byte  or  16  bit 
format)  except  for  the  COUNT  and  DELTA  data  sets  which  are 
1*4  (cell  values  are  represented  as  an  integer  value  stored  in  a 
four  byte  or  32-bit  format).  The  maximum  integer  value  that  can 
be  stored  in  a  two  byte  format  is  32767,  so  the  highest  integer 
value  of  a  cell  in  a  map  stored  in  this  fashion  is  32767.  The 
maximum  integer  value  that  can  be  stored  in  a  four  byte  format 
is  2147483647,  and  therefore  a  much  larger  integer  cell  value  is 
possible.  The  COUNT  and  DELTA  data  sets  often  generate 
integer  cell  values  much  higher  than  32767  and  therefore 
require  the  1*4  format.  This  becomes  relevant  when  discussing 
the  interface  with  MOSS/MAPS. 

Restart  capability 

The  COUNT  and  DIRECT  programs  contain  a  restart  capability 
that  permits  resumption  of  processing  if  the  program  stopped 
because  of  insufficient  disk  space,  or  any  other  reason.  Enter  a 
1  (one)  to  restart  the  process.  Enter  a  0  (zero)  to  begin  the 
program  afresh.  This  capability  has  not  been  tested. 
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Information  messages 

Many  of  the  programs  are  iterative  and  report  the  completion  of 
each  iteration  for  user  information.  CONCAVE  prints  a  message 
as  it  initializes  each  box,  and  SEED  reports  how  many  basin 
starts  it  seeded.  In  addition  various  error  messages  appear. 

Program  Limitations 

The  programs  limit  the  number  of  samples  (columns)  to  4000. 
CONCAVE  limits  BOXNL  and  BOXNS  to  400,  or  the  row  and 
column  size  of  a  rowing  window.  The  PPT  ABLE  program  limits 
its  pour  point  table  entries  to  2000.  If  you  encounter  the  pour 
point  table  limit  indicated  by  a  "STOP  72"  message,  it  is  pos- 
sible to  process  the  data  set  in  pieces.  The  RSNAP  program 
limits  the  threshold  distance  to  50  cells  in  the  x  and  y  direction 
and  does  not  permit  movement  of  more  than  999  points. 

IVIOSS/MAPS  Interfacing 

The  Hydrologic  Terrain  Analysis  Software  (HTAS)  interfaces 
with  MOSS/MAPS  as  follows. 

I     The  MAPSTOIDIMS  command  permits  conversion  of  a 
MAPS  map  into  IDIMS  format  (a  sequential  access 
format).  A  CONVERT  program  converts  the  IDIMS  data 
into  the  direct  access  format  used  by  the  HTAS. 

I     Use  a  reverse  path  to  bring  an  HTAS  processed  map  into 
MAPS.  The  CONVERT  program  transfers  the  MAPS  file  it 
from  "direct"  to  IDIMS  format,  and  then  the  MAPS  IMPORT 
command,  WORD  option,  imports  the  map  into 
MOSS/MAPS. 


NOTE: 

All  the  raster  maps  produced  by  the  HTAS  use  the  1*2  (integer  *  2 
bytes)  format,  with  the  exception  of  the  flow  accumulation  map 
which  has  the  1*4  (integer  *  4  bytes)  format.  Currently,  the  MAPS 
IMPORT  command  only  imports  1*2  data,  with  the  highest 
allowable  integer  number  being  32768.  There  is  currently  an 
improved  version  of  IMPORT,  included  with  the  HTAS  package, 
which  permits  conversion  of  1*4  files  into  MAPS  format  as  real 
numbers  rather  than  integer  cell  value  maps. 
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Automated  Hydrologic  Terrain 

Analysis  Concept 


Influence  ef  terrain  morphology  on  water  flow 

As  rainwater  falls  on  the  landscape,  the  morphology  of  the 
terrain  is  the  major,  if  not  the  only  determinant  of  which  way 
the  water  flows.  For  the  purpose  of  this  discussion,  landscape 
morphology  includes  all  the  macroscopic  details  of  the  shape 
of  the  surface.  Surface  factors  that  affect  water  flow  also  include 
the  obstructive  effect  of  vegetation  cover,  the  infiltrative  capac- 
ity and  moisture  content  of  soil,  the  presence  of  cracks  or 
fissures  in  rocks  and  artificially  made  obstacles  and  barriers. 
Although  some  of  their  aspects  could  be  related  to  surface 
morphology,  these  factors  include  physical,  biological  and 
chemical  processes  which,  while  influencing  water  flow,  are 
not  directly  related  to  the  gross  shape  of  the  terrain.  This 
manual  does  not  discuss  these  factors  or  the  factors  affecting 
water  flow  that  relate  to  the  processes  occurring  underneath 
the  surface. 

Elevation  data  and  their  resolution 

Elevation  contour  maps,  digital  elevation  models  (DEMs)  and 
digital  terrain  models  of  various  scales  and  resolutions,  provide 
basic  information  about  the  general,  macroscopic  trends  in  the 
shape  of  the  landscape  surfaces.  A  typical  USGS  digital  eleva- 
tion model  has  a  resolution  of  30x30  meters  and  clearly  loses  a 
significant  amount  of  information  about  the  detail  of  the  surface 
when  only  one  number  represents  a  square  area  with  a  side  of 
98.4  feet.  Therefore,  all  information  obtained  through  HTAS 
processing  with  data  of  this  resolution  provides  only  general 
trends  and  approximations  of  what  actually  occurs,  or  what 
might  occur.  The  HTAS  only  provides  information  on  water 
flow  based  on  the  gross  shape  of  the  landscape,  within  the 
resolution  of  the  data  provided,  and  with  the  assumption  that 
an  unlimited  supply  of  water  runs  over  an  essentially  imperme- 
able terrain.  To  perform  realistic  modeling  of  surface  hydrologic 
processes  we  must  combine  the  information  obtained  using 
HTAS  with  other  data  and  processing. 
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Filling  depressions  on  the  landscape 

The  landscape  often  has  depressions  where  the  water  pools 
until  it  finally  overflows  and  then  continues  to  pour  in  the 
direction  determined  by  the  terrain  and  other  landscape  factors. 
Imagine  a  tilted  cup  into  which  you  are  pouring  water.  If  we 
have  an  unlimited  supply  of  water,  the  cup  eventually  fills  and 
finally  water  flows  out  at  the  lowest  point  on  the  edge  of  the 
cup.  Therefore,  to  route  water  flow  beyond  depressions  using 
the  HTAS  it  is  first  necessary  to  fill  them. 


NOTE: 

In  reality,  such  complete  filling  of  a  sink  might  never  occur.  In 
dry  areas  with  slight  rainfall  the  water  might  not  overflow  a 
large  depression  but  pond  in  the  depression  and  infiltrate  the 
soil.  Such  depressions,  depending  on  the  local  soil  and  geology, 
might  become  recharge  locations  for  an  aquifer  or  sediment 
collection  areas. 


The  HTAS  contains  programs  which  fill  depressions  on  the 
surface  of  a  DEM.  Depressions  found  in  the  DEMs  reflect 
natural  conditions  on  the  ground  or  imperfections  in  the  model 
itself.  Generally,  sinks  or  pits  where  only  one  cell  has  a  value 
lower  than  the  surrounding  cells  reflects  an  imperfection  in  the 
DEM,  however,  there  are  cases,  for  example  in  glacial  or  karst 
areas,  where  these  pits  may  reflect  natural  features.  Multi-cell 
depressions  are  more  likely  to  reflect  natural  features.  First  hand 
knowledge  of  field  conditions  is  very  important  when  interpret- 
ing the  results  of  the  depression-filling  process. 

The  program  that  deals  with  single  cell  sinks  is  FILLSNGL, 
which  raises  the  value  of  the  depressed  cell  to  the  value  of  its 
lowest  immediate  neighbor  in  a  3x3  cell  "window".  The  de- 
pressed cell  is  at  the  center  of  the  window  which  includes  its 
nine  immediate  neighbors.  Figure  1  shows  single  cell  depres- 
sions that  were  found  and  filled  using  FILLSNGL  in  the  Digital 
Elevation  Model  for  the  7.5  minute  Sagebrush  Hill  quadrangle. 
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Figure  1:  Stream  network  for  a  section  of  the  Sagebrush  Hill  Quadrangle,  Colorado,  with 
single  cell  sinks  filled  using  the  FILLSNGL  program. 
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The  second  program  used  to  fill  depressions,  CONCAVE, 
deals  with  large,  multi-celled  depressions  in  the  model.  The 
program  requires  the  user  to  specify  the  size  of  the  window  and 
an  overlap  area  in  terms  of  rows  and  columns.  The  size  of  the 
depression  must  be  smaller  than  the  size  of  the  overlap  for 
program  to  fill  the  depression.  Figure  2  shows  larger  depressions 
found  and  filled  using  CONCAVE  for  the  same  quadrangle. 


Figure  2:  Stream  network  for  a  section  of  the  Sagebrush  Hill  Quadrangle,  Colorado, 
ivith  large  depressions  filled  using  the  CONCAVE  program. 
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Determining  flow  direction 

Once  HTAS  programs  fill  the  depressions  in  the  DEM,  it  is 
possible  to  route  the  water  completely  across  the  landscape  and 
discover  the  direction  in  which  it  flows  on  the  depressionless 
surface.  The  program  DIRECT  provided  with  the  HTAS  calcu- 
lates flow  directions  on  the  depressionless  surface  while  paying 
special  attention  to  connectivity  between  cells.  One  of  the  eight 
possible  directions  to  adjacent  cells  at  the  north,  northeast,  east, 
southeast,  south,  southwest,  west,  and  northwest  of  a  particular 
cell  is  assigned  to  that  cell.  Based  on  this  information  it  be- 
comes possible  to  calculate  flow  accumulation. 

Flow  accumulation  and  drainage  pattern  development 

Drainage  patterns  develop  when  the  rate  of  rainfall  exceeds 
the  rate  of  infiltration,  causing  the  depressions  on  the  surface 
to  fill  and  overflow  forming  a  thin  sheet  of  water  that  flows 
across  the  landscape.  This  runoff  takes  the  path  of  least  resis- 
tance and  gains  velocity.   It  picks  up  soil  particles  and  creates 
small  channels,  so  that  the  surface  of  the  ground  exhibits,  on  a 
miniature  scale,  the  drainage  pattern  of  a  watershed  by  forming 
rills.  Finally  a  fluctuating  drainage  pattern  develops  throughout 
the  entire  catchment  which  includes  rills,  channels,  ephemeral 
and  perennial  streams,  and  finally  rivers. 

The  flow  accumulation  program  COUNT  in  the  HTAS  permits 
approximate  delineation  of  drainage  patterns  from  DEMs  based 
on  flow  directions  established  by  the  information  on  the  gross 
shape  of  the  landscape.  Using  the  flow  direction  data  set  the 
program  evaluates  how  many  other  cells  are  connected  and 
flow  into  a  particular  cell.  Therefore,  the  value  of  each  cell  in  a 
flow  accumulation,  or  flow  count  map,  is  equivalent  to  the 
number  of  cells  that  flow  into  this  particular  cell. 

Since  a  DEM  provides  only  a  gross  representation  of  the  land- 
scape, HTAS  identifies  flow  patterns  based  only  on  the  general 
shape  of  the  terrain.  In  addition,  it  does  not  consider  the  influ- 
ence of  other  environmental  factors,  including  microrelief. 
The  program,  however,  predicts  flow  patterns  with  surprising 
accuracy,  particularly  in  areas  of  high  relief.  On  large  flat  areas 
the  influence  of  minor  variations  in  the  shape  of  the  land  which 
do  not  show  up  in  DEMs  becomes  more  important,  and  the 
program  is  subject  to  a  greater  degree  of  error. 

Figure  3  shows  the  drainage  pattern  detected  by  the  program  in 
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a  high  relief  terrain  in  a  section  of  the  Sagebrush  Hill  quad- 
rangle. When  the  stream  network  information  derived  from  a 
USGS  Digital  Line  Graph  overlays  the  modeled  patterns,  we  see 
the  close  correspondence  between  the  two  data  sets. 


Figure  3:  DLG  stream  data  (dark  continuous  lines)  overlay  the  modeled  drainage  patterns 
which  include  potential  and  actual  channels  from  the  Sagebrush  Hill  quadrangle. 

The  modeled  data  show  locations  of  a  likely  drainage  pattern 
based  on  gross  morphology  of  the  surface,  and  should  be 
treated  with  caution  when  the  areas  modeled  have  low  relief. 
The  modeled  data  do  not  replace  stream  network  information 
obtained  from  aerial  photographic  surveys,  but  the  models  are 
very  useful  for  areas  where  there  are  no  other  data  sources 
available. 
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Watershed  boundary  delineations 

A  watershed,  drainage  basin,  catchment,  or  a  contributing  area 
is  an  area  of  land  where  water  flows  toward  a  particular  loca- 
tion. That  location  can  be  a  dam,  a  lake,  or  simply  a  pour  point 
for  a  particular  watershed.   You  can  divide  a  watershed  into 
subwatersheds.  Generally,  this  division  exists  at  a  point  where 
a  tributary  stream  enters  another  stream  in  a  drainage  network. 
The  total  of  specific  subwatersheds  for  all  of  the  tributaries  of  a 
particular  main  stream  constitutes  the  drainage  basin  for  the 
main  stream.  It  is,  however,  possible  to  delineate  a  catchment 
area  for  any  point  on  the  landscape,  as  long  as  the  shape  of  the 
land  is  such  that  there  is  water  flowing  towards  that  point.  The 
capabilities  of  HTAS  do  not  limit  analysis  to  pour  points  at  the 
nodes  of  a  drainage  network. 

Watershed  delineation  with  the  HTAS  identifies  all  the  cells 
flowing  into  a  particular  pour  point,  line  (e.g.,  dam),  or  an  area 
(e.g.,  reservoir)  as  belonging  to  a  single  watershed.  Based  on 
the  filled  DEM,  flow  direction,  and  flow  accumulation  data  sets 
obtained  through  prior  processing,the  HTA  Software  permits 
delineation  of  watershed  and  subwatershed  boundaries  both  in 
a  manual  and  automated  fashion.  The  automated  method 
selects  pour  points  for  possible  watersheds  based  on  a  user- 
selected  value  of  flow  accumulation,  which  is  called  a  threshold 
value.  The  manual  method  permits  the  user  to  select  pour 
points,  linear  or  polygonal  (area)  features,  for  which  they  desire 
to  delineate  a  watershed.  Both  cases  require  creation  of  a 
"seed"  or  "starter"  map;  the  word  "seed"  here  being  synony- 
mous with  the  representation  in  a  cell  map  of  a  point,  line  or 
polygon  for  which  a  watershed  is  delineated. 
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Figures  4,  5,  6  show  watersheds  generated  for  points,  lines  and 
polygons  representing  watershed  outlets,  hypotlietical  dams 
and  reservoirs.  Figure  7  shows  the  watersheds  generated  for  a 
threshold  value  of  3000  for  the  Sagebrush  Hill  quadrangle.  The 
delineated  watersheds  do  not  reach  the  edges  of  the  DEM,  an 
option  offered  with  the  HTAS. 


Figure  4.  Watersheds  delineated  for  user-selected  pour  points. 
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Figure  5.  Watersheds  delineated  for  user-selected  linear  features. 
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Figure  6.  Watersheds  delineated  for  user-selected  reservoirs. 
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Figure  7.  Automatically  delineated  subwatersheds  using  the  threshold  values  of  3000. 
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Overland  flow  path  delineation 

Another  option  of  the  HTAS  calculates  and  delineates  the 
steepest  overland  path  from  a  point  or  a  set  of  points  on  the 
landscape.  The  program  OVERLAND  utilizes  a  flow  direction 
file  and  a  mask  file  similar  to  a  seed  file,  which  identifies  points 
on  the  landscape  from  which  to  delineate  the  path.  Figure  8 
provides  an  example  of  overland  paths  calculated  for  three 
points  in  the  Sagebrush  Hill  quadrangle. 

Mapping  local  minima  and  another  method  for  development 
of  a  drainage  network 

The  program  MINIM,  when  used  on  the  elevation  model, 
permits  identification  of  the  local  minima,  or  cells  that  have  a 
lower  elevation  value  than  the  surrounding  cells  in  a  3x3  cell 
window.  This  program  also  permits  identification  of  flat  areas 
where  the  central  cell  has  the  same  value  as  all  the  surrounding 
cells  in  a  3x3  window.  When  the  data  set  produced  by  MINIM 
is  processed  using  OVERLAND,  the  OVERLAND  program  traces 
the  overiand  paths  for  all  these  cells  using  the  flow  direction 
data  set.  The  program  marks  all  the  cells  into  which  a  cell  of 
the  MINIM  set  drains  until  it  reaches  another  cell  of  the  MINIM 
set  that  is  on  the  path.  The  result  is  the  image  of  a  drainage 
network  based  on  the  elevation  data  set.  This  is  a  second 
method  for  creating  a  drainage  network,  since  the  drainage 
pattern  can  also  be  delineated  using  the  flow  accumulation  data 
set  as  described  before. 
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Figure  8.  Steepest  overland  paths  from  point  sources  calculated  using  the  OVERLAND 
program  join  with  streams. 
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Figure  9.  The  drainage  network  data  set  generated  using  the  MINIM  and  OVERLAND  pro- 
grams with  the  DIG  stream  data  OVERLAND. 
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I.  Convert  and  transfer  files  between 

MOSS/MAPS  and  the  H1AS 


GIS  Vector 
File 


MAPS  RASTERIZE 

or 

MOSS  POLYCELL 


The  conversion  and  transfer  of  files  between  GIS  and  HTAS 
requires  several  steps  and  is  somewhat  cumbersome.  In  the 
future,  HTAS  will  potentially  integrate  MOSS/MAPS  capabilities 
which  negates  the  need  for  file  transfer  between  the  two 
formats.  Currently  conversion  and  transfer  use  MAPSTOIDIMS 
and  CONVERT  commands,  utilities  which  are  executed  at  the 
Primos  level,  and  the  MAPS  IMPORT  command. 

Maps  exist  or  can  be  generated  in  MOSS  vector  format  and 
MAPS  raster  format.  Use  the  POLYCELL  command  in  MOSS  or 
the  RASTERIZE  command  in  MAPS  to  convert  maps  from  MOSS 
vector  format  to  MAPS  raster  format.  The  HTAS  works  only  with 
raster  data.  The  MOSS  User's  Manual  discusses  the  POLYCELL 
and  RASTERIZE  commands  in  detail;  please  refer  to  that  docu- 
ment for  additional  information. 


GIS  Raster 
File 


CIS  Raster 
File 


MAPSTOIDIMS 


MAPS  IMPORT 


Sequential 
Access  File 


CONVERT 


Sequential 
Access  File 


CONVERT 


Direct  Access 
File 


Direct  Access 
File 


f  sw*\\iW!,  -*^B 
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STEP  1      Converting  MAPS  format  maps  into  sequentiai  access 
raster  format  (IDiMS)  maps. 

Use  the  MAPSTOIDIMS  program  to  convert  a  MAPS  map  into  a 
sequential  access  raster  transfer  format  (IDIMS  format  named 
for  the  image  processing  system,  IDIMS,  which  BLM  has  been 
using).  MAPSTOIDIMS,  an  operating  system  (PRIMOS)  level 
utility,  converts  raster  files  in  MAPS  format  into  the  sequential 
access  raster  transfer  format  (IDIMS  format).  Originally  written 
to  convert  and  write  files  to  tape,  it  also  converts  and  writes 
files  to  disk  by  using  defaults.  MAPSTOIDIMS  creates  files  with 
.16  and  .TH  suffixes  after  the  name  of  the  input  file.  The  .16 
files  contain  the  raster  files  used  for  processing,  while  the  .TH 
files  contain  information  about  various  parameters  of  the  master 
file  such  as  the  number  of  rows  and  columns  or  coordinates  of 
the  Minimum  Bounding  Rectangle  (MBR)  of  the  map  (the 
window  of  the  map)  to  which  the  user  may  refer. 


WARNING: 

MAPSTOIDIMS  works  only  with  l6-bit  (word  length)  integer 
data;  however,  some  of  the  MAPS  real  number  files  might  be  32- 
bit  (double  word)  data.  To  correct  this  problem  use  the  FUNC- 
TION command  INTEGER  option  in  MAPS  to  convert  the  real 
number  32-bit  map  into  an  integer  1 6-bit  map. 
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OK,  HflPSTOIDIMS 

MAPS  TO  ID  IMS  CONUERSION 

Enter  macro  file  name  to  execute  file  copy 

<  Default  is  TMCOPV  ) 
<CR> 

Enter  tope  drive  name  MTO  or  MT1  (  Default  =  MTO  ) 
<CR> 

Enter  tape  drive  density  (800,  1600  or  6250;  Default  =  1600  ) 

Houi  many  maps  are  already  on  this  tope  ? 
:<CR> 

Please  enter  the  DIRECTORV  (under  MOSSDflTfl)  or  NORK 

and  the  .DT  (project  master  or  POLVGON  ; 

i.e.,  UOLSflG  UOLSRG.DT  or  WORK  POLYGON. DT 
:  UOLSRO  32  U0LSRG32.DT 
Enter  map  name  (  CR  =  End  ) 
:  DEMFILE  (name  of  a  tIFIPS  map) 

Map-  DEMFILE       has  been  processed. 
Enter  map  name  (  CR  =  End  ) 
:<CR> 
*  INSTRUCT  I OMS  TO  COMPLETE  MflPSTOIDIMS  * 

1)  List  out  the  .TH  file  created  in  the 
same  directory  that  the  MflPSTOIDIMS  was 
executed  (SPOOL  or  SLIST  on  the  PRIME). 
This  file  contains  the  MBR,  cell  size, 
and  projection  criteria  of  the  file, 
and  needs  to  accompany  the  TAPE 
transferred  to  I  DIMS 

2)  Run  the  macro  called  TMCOPV  to  transfer 
the  cell  file  to  magnetic  tape: 

t      Mount  the  tape 

•       Run  TMCOPV  (CPL  TMCOPV) 

MflPSTOIDIMS  completed. 
****  5J0P 


NOTE: 

The  necessity  for  conversion  and  transfer  of  map  data  via  the 
sequential  access  raster  transfer  (IDIMS)  format  rests  on  the  fact 
that  at  this  stage  the  system  is  available  only  on  PRIMOS,  and  no 
programming  additions,  except  for  minor  changes  in  the  MAPS 
IMPORT  command  are  anticipated  at  this  time.  In  the  future, 
HTAS  will  contain  more  direct  interfaces  should  funds  become 
available  for  further  improvements. 
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STEP  2      Converting  sequential  access  raster  transfer  (IDIiVIS) 
(ormat  maps  into  direct  access  format  maps  used  in  tlie 
HTAS. 

After  converting  the  data  into  the  IDIMS  format,  it  must  be 
converted  again  into  direct  access  format  which  the  HTAS  uses. 
Use  the  CONVERT  program  as  follows  to  accomplish  this  task. 


OK,  R  CONUERT 

Conversion  desired  (0=1*2  transfer- to-direct,    1=1*2  or 
1*4  direct-to-transfer): 
0 

Transfer  DEM  input  image  name: 
DEH.ttflP.16 

Direct  access  output  DEM  image: 
DEM. MRP 

Lines  in  image: 
472 

Samples  in  image: 

36g 

Converted        200  I ines 
Converted        400  I Ines 
4t)ic)ii)|[  STOP 


STEP  3      Processing  direct  access  format  maps  in  tlie  HTAS. 

Once  the  data  are  converted  into  the  direct  access  format,  any  of 
the  HTAS  programs  can  use  the  data  set.  The  HTAS  Reference 
Manual  contains  an  alphabetic  listing  of  each  program  and  a 
description  of  specific  hydrologic  modeling  applications  appear 
in  the  following  parts  of  the  HTAS  User's  Manual. 


STEP  4      Converting  direct  access  format  maps  produced  by  tiie 
HTAS  processing  into  sequentiai  access  raster  transfer 
(IDiiVIS)  format. 

To  bring  the  data  back  into  MAPS,  convert  the  map  into  IDIMS 
format  using  the  program  CONVERT  again.  An  example  of 
converting  from  direct  to  sequential  access  raster  transfer 
(IDIMS)  format  follows. 
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OK,   R  CONUERT 

Conversion  desired  (0=1*2  transfer-to-direct,  1=1*2  or 
1*4  direct-to-transfer): 
1 

Direct  access  input  image  name: 
SHED.nRP 

Transfer  output  image  name: 
SHED.nRP. 16 

Direct- to- transfer  data  type  (2=1*2,  4=1*4): 
2 

Lines  in  image: 
472 

Samples  In  image: 
360 

Converted        200  I ines 

Converted        400  I ines 
****  STOP 


NOTE: 

The  flow  accumulation  set  produced  by  the  COUNT  program 
produces  an  1*4  map  file.  To  convert  1*4  files  into  the  sequential 
access  raster  transfer  (IDIMS)  format,  select  option  4  in  the  "Direct- 
to-transfer"  data  type  prompt. 


STEP  5      Importing  sequential  access  raster  transfer  (IDIMS)  format 
maps  Into  the  Map  Analysis  and  Processing  System  (MAPS) 
format. 

After  the  map  file  is  converted  to  the  IDIMS  format,  use  the  MAPS 
IMPORT  command  to  bring  the  map  into  the  MAPS  GIS.  The 
original  IMPORT  permits  importation  of  word  length  or  1*2  maps 
only.  The  new  version  included  with  this  software  package 
contains  the  capability  to  bring  "double  word"  or  1*4  maps  as  real 
number  maps  into  the  MAPS  GIS. 

Below  is  an  example  of  how  the  MAPS  IMPORT  command 
works  with  the  1*2  maps;  for  the  import  1*4  maps,  use  the 
"DOUBLEWORD"  instead  of  the  "WORD"  option  in  the  new 
version  of  IMPORT.  The  "DOUBLEWORD"  option  permits  impor- 
tation of  type  8  maps  only. 


User's  Manual  Part  1 


3D1HTAS  User  8  Manual 


IMPORT  UflTERSHED.16  FORMflT    IDIMS  UORD  FOR  HRTERSHED  TVPE  7 

ENTER  THE  FOLLOWING   INFORMflTION  FOR   INPUT  MRP 

NUMBER  OF  ROUS  ? 
472 

NUMBER  OF  COLUMNS  ? 
360 

CELL  HEIGHT  (IN  METERS)  ? 
30 

CELL  WIDTH  (IN  METERS)  ? 
30 

ENTER  THE  MINIMUM  BOUNDING  RECTANGLE 

MINIMUM  X  (WEST)  ? 
702750 

MAXIMUM  X  (ERST)  ? 
713820 

MINIMUM  V  (SOUTH)  ? 
4416390 

MRXIMUM  V  (NORTH)  ? 
4430550 

ENTER  MRP  PROJECTION   INFORMATION 

PROJECT  1 0N(0-20)  ? 
1 

ELLIPS0ID(0-19)  ? 
0 

LONGITUDE  OF  RNV  POINT  WITHIN  THE  UTM  ZONE  ? 
-109 

LATITUDE  OF  ANV  POINT  WITHIN  UTM  ZONE       12  ? 
35 

ENTER  MRP  DESCRIPTION   INFORMflTION 

MAP  DESCRIPTION  ? 
URTERSHEDS  BASED  ON  THRESHOLD  OF  3000 

STUDV  RREA  ? 
UOLSRO 

PROJECTION  DESCRIPTION  ? 

uin 

DRTE  ? 
10/01/92 

SOURCE  ? 
J.BLflSZCZVNSKI 

UINTRGE  ? 
1992 
OK     IMPORTED  FOR  URTERSHED 

NOTE: 

Once  the  map  is  imported  iJito  the  MAPS  CIS  it  becomes  available 
for  use  as  a  data  layer  in  analysis  and  cartographic  modeling 
performed  using  MAPS  capabilities.  The  final  chapter  of  this 
manual  briefly  discusses  these  potential  applications. 
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Terrain  Data 

Conditioning  Procedures 


S.  Jenson  refers  to  the  aspects  of  HTAS  processing  that  result  in 
a  depressionless  DEM,  a  flow  direction  data  set  and  a  flow 
accumulation  data  set,  as  the  conditioning  procedures  because 
they  produce  the  data  sets  that  are  of  general  utility  for  all 
subsequent  modeling. 


STEP  1      Developing  a  depressioniess  DEM 


DEM 


HLLSNGL 


DEM 
with  Single  Cell 
Depressions  Filled 


CONCAVE 


DEM 
with  Multi-Cell 
Depressions  Filled 


DIRECT 


The  Flow  Direction 
Data  Set 


The  Flow  Accumulation 
Data  Set 


HTAS  provides  two  programs  which  fill  depressions  in  a 
DEM  to  determine  the  flow  of  water  to  the  edge  of  the  data  set 
across  the  modeled  terrain.  The  first  program,  FILLSNGL  fills 
only  the  single  celled  depressions  which  are  usually  a  result  of 
errors  in  the  DEM.  The  second  program,  CONCAVE,  fills  multi- 
cell  depressions  that  usually  represent  the  actual  shape  of 
the  terrain. 

FILLSNGL  processes  data  using  a  3x3  cell  window,  and  raises 
the  value  of  the  center  cell  to  the  value  of  the  lowest  of  its 
neighbors.  CONCAVE  processes  the  data  using  overlapping 
windows  of  user  specified  dimensions,  i.e.,  size  of  the  window 
in  lines  (rows)  and  samples  (columns)  and  size  of  the  overlap. 
The  size  of  the  overlap  must  be  greater  than  that  of  the  depres- 
sion being  filled. 

This  process  results  in  a  depressionless  DEM  where  each  cell  is 
a  part  of  at  least  one  monotonically  decreasing  path  of  cells  that 
leads  to  the  edge  of  the  data  set.  This  path  consists  of  cells  that 
are  adjacent  horizontally,  vertically  or  diagonally  in  an  eight- 
way  connectedness  and  that  steadily  decrease  in  elevation 
value  toward  the  edge  of  the  data  set. 

Since  both  FILLSNGL  and  CONCAVE  read  from  and  write  to 
the  same  file,  keep  a  copy  of  the  original  DEM  in  the  special 
case  where  flow  routing  with  depressions  present  becomes 
important.  While  CONCAVE  fills  all  the  depressions,  it  is  useful 
to  execute  FILLSNGL  first  and  keep  a  copy  of  the  product  to 
clearly  distinguishing  between  single-cell  and  multi-cell 
depressions. 

For  more  details  on  the  programs  refer  to  the  HTAS  User's 
Reference  Manual. 
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STEP  1A    Executing  FILLSNGL 


To  execute  FILLSNGL  start  at  that  directory  where  the  execut- 
able version  of  the  program  FILLSNGL.RUN  resides  by  typing  R 
FILLSNGL,  or  identify  the  path  to  the  program  from  the  current 
directory  by  typing  R  (pathname)>FILLSNGL.  The  inputs  to  the 
program  include  the  number  of  lines  or  rows  (NL),  the  number 
of  samples  or  columns  (NS),  the  name  of  the  unprocessed, 
original  elevation  file  in  single  quotes  (TOPO  FILE)  which  is 
the  input  and  the  file  to  which  the  output  is  written,  and  the 
name  of  the  system  in  single  quotes,  i.e.,  'PRIME'  (SYSTEM). 


OK,  R  FILLSNGL 

ENTER  NL,NS,LEU  FILE,SVSTEM(UNIX,UMS,PRiriE) 
472,365, 'DEn.nflP', 'PRIME' 
Open  as  a  Prime  direct  access  file,  ns  = 
DEM. MAP 

opened  successfully 
)iitt*  STOP 


369 


^  \;;^>^^'"\'^A^5?'^^*"'''^V^'*^'^'^^^ 


NOTE: 

FILLSNGL  reads  from  and  writes  to  the  same  file,  therefore  keep 
the  original  file  in  the  HTAS  format  and  use  a  copy  for  FILLSNGL 
processing. 
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STEP  IB   Executing  CONCAVE 

To  execute  CONCAVE  start  at  that  directory  where  the  execut- 
able version  of  the  program  CONCAVE. R  resides  by  typing  R 
CONCAVE,  or  identify  the  path  to  the  program  from  the  current 
directory  by  typing  R  (pathname)>CONCAVE.  The  inputs  to  the 
program  include  the  number  of  lines  or  rows  (NL),  the  number 
of  samples  or  columns  (NS),  window  size  in  lines  or  rows 
(BOXNL),  window  size  in  samples  or  columns  (BOXNS), 
overiap  between  windows  (OVERLAP),  the  name  of  the  eleva- 
tion file  in  single  quotes  (LEV  FILE),  and  the  name  of  the 
system  in  single  quotes,  i.e.,  'PRIME'  (SYSTEM). 


NOTE: 

CONCAVE  reads  from  and  writes  to  the  same  file,  therefore  keep 
the  original  file  in  the  HTAS  format  and  use  a  copy  for  CON- 
CA  VE  processing. 


OK,  R  CONCRUE 

ENTER  NL,NS, BOXNL, BOXNS,  OUERLflP, LEU  FILE,SVSTEf1<UNIX,UNS,PRmE) 
472,36Q,  '100, 100,60, 'CONCRUE. HRP', 'PRIHE' 
System  =  PRIME 

Open  as  a  Prime  direct  access  file,  ns  =     369 
SGML. CONCRUE 
opened  successful ly 


SL  = 

1 

SL  = 

41 

SL  = 

81 

SL  = 

121 

SL  = 

161 

SL  = 

201 

SL  = 

241 

SL  = 

281 

SL  = 

321 

SL  = 

361 

SL  = 

401 

SL  = 

441 
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STEP  2      Developing  a  flow  direction  data  set. 

The  HTAS  program  DIRECT  creates  the  flow  direction  data  set 
and  encodes,  for  each  cell  of  the  DEM,  the  direction  in  which 
the  water  flows  out  of  that  cell.  The  following  schematic  illus- 
trates the  directions  of  water  flow  as  they  are  encoded  in  the 
flow  direction  map. 


64(NW) 

128(N) 

1(NE) 

32(W) 

center 

2(E) 

16(SW) 

8(S) 

4(SE) 

Although  the  flow  direction  data  set  is  created  as  part  of  the 
conditioning  stage  for  delineating  watersheds,  its  also  assists  in 
visualizing  general  flow  patterns  on  the  landscape  when  dis- 
played. Differing  from  the  maps  created  by  the  ASPECT  com- 
mand in  MAPS,  this  data  set  describes  the  direction  of  water 
flow  in  the  eight-connected  cells  within  a  nine  cell  neighbor- 
hood, rather  than  calculating  the  directions  in  which  the  slope 
of  the  cell  faces.  Apply  the  DIRECT  command  to  either  an 
unfilled  DEM  or  to  a  filled  DEM  (see  FILLSNGL  and  CONCAVE). 
When  applied  to  an  unfilled  or  partially  filled  DEM,  the  flow 
direction  data  set  contains  negative  numbers  to  represent  the 
cells  for  which  the  flow  could  not  be  calculated.  When  applied 
to  a  completely  filled  DEM,  the  program  assigns  positive  flow 
directions  to  each  cell,  which  in  turn  allows  calculation  of  a 
monotonically  decreasing  continuous  path  from  each  cell  to  the 
edge  of  the  DEM.  For  more  details  on  processing  in  DIRECT 
see  the  appropriate  section  in  the  HTAS  User's  Manual. 
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To  execute  DIRECT  start  at  that  directory  where  the  executable 
version  of  the  program   DIRECT.RUN  resides  by  typing  R 
DIRECT,  or  identify  the  path  to  the  program  from  the  current 
directory  by  typing  R  (pathname)>DIRECT.  The  inputs  to  the 
program  include  the  number  of  lines  or  rows  (NL),  the  number 
of  samples  or  columns  (NS),  a  restart  option  if  processing 
inten-upted  (CODE  1  TO  RESTART),  the  name  of  the  elevation 
file  in  single  quotes  (LEV  FILE),  the  name  of  the  flow  direction 
file  in  single  quotes  (DIR  FILE),  and  the  name  of  the  system  in 
single  quotes,  i.e.,  'PRIME'  (SYSTEM). 

OK,   R  DIRECT 
ENTER  ML, NS, CODE  1  TO  RESTART, LEU  FILE,DIR  FILE,SVSTEn<U 

NIX, UnS, PRIME) 
472,36Q,0, 'DEM. MRP', 'DIRECT. HRP', 'PRIME' 
Open  ELEU  and  DIR  files,   ns  =  369 

DEM. MRP 

opened  successful ly 
DIRECT. MRP 

opened  successful ly 
DOUMURRD  PASS    1FIRSTL  2LRSTL    471 
UPWRRD  PRSS        2FIRSTL  7LRSTL    465 
DOWNWARD  PRSS    3FIRSTL  7LRSTL    461 
COULD  NOT  SOLUE  FOR  RLL  CELLS 
****  STOP 


STEP  3      Developing  a  flow  accumulation  data  set. 

Create  the  flow  accumulation  data  set  by  executing  the  program 
COUNT.   COUNT  computes  flow  accumulation  values  from 
flow  direction  values.  Each  cell  in  the  COUNT  data  set  receives 
a  value  equal  to  the  number  of  cells  that  flow  into  it.  Cells 
having  the  flow  accumulation  value  of  zero  (to  which  no  other 
cells  flow)  generally  correspond  to  the  pattern  of  ridges. 

Because  all  cells  in  a  depressionless  DEM  (a  product  of 
FILLSNGL  and/or  CONCAVE)  have  a  path  to  the  data  set  edge, 
the  pattern  formed  by  highlighting  cells  with  values  higher  than 
a  given  threshold  value  delineates  a  fully  connected  drainage 
network.  As  the  threshold  value  is  decreased,  the  density  of  the 
drainage  network  increases.  COUNT  is  a  part  of  the  condition- 
ing stage  for  delineating  watersheds,  however  its  product  also 
has  a  value  for  delineating  drainage  networks  and  potential 
stream  flow  modeling. 
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To  execute  COUNT  start  at  the  directory  where  the  executable 
version  of  the  program  COUNT.R  resides  by  typing  R  COUNT, 
or  identify  the  path  to  the  program  from  the  current  directory 
by  typing  R  (pathname)>COUNT.  The  inputs  to  the  program 
include  the  number  of  lines  or  rows  (NL),  the  number  of 
samples  or  columns  (NS),  a  restart  option  if  processing  inter- 
rupted (CODE  1  TO  RESTART),  the  name  of  flow  direction  file 
in  single  quotes  (DIR  FILE),  the  name  of  the  flow  accumulation 
file  in  single  quotes  (COUNT  FILE),  and  the  name  of  the  system 
in  single  quotes,  i.e.,  'PRIME'  (SYSTEM). 


OK,  R  COUNT 

ENTER  NL,NS,CODE   1  TO  RESTART, DIR  FILE, COUNT 

FILE,SVSTEN(UNIX,UMS,PRmE) 
472,360,0, 'DIRECT. nflP', 'COUNT. MRP', 'PRIME' 


DOUNURRD  PRSS 
UPWARD  PASS 
DOWNWARD  PASS 
UPWARD  PASS 
DOWNWARD  PRSS 
UPWARD  PRSS 
DOWNWARD  PASS 
UPWARD  PASS 
DOWNWARD  PASS 
UPWARD  PASS 
DOWNWARD  PASS 
UPWARD  PRSS 
DOWNWARD  PASS 
UPWARD  PASS 

****    5JQP 
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12 


14 


1 


11 


13 
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Watersheds: 

Manual  Delineation 


Display  of  Flow 
J  Accumulation  Drainages  j 
as  Backdrop  to  Guide 
Placement  of  Pour 
iLocations  for  Watersheds^ 


MOSS  GENERATE 


Vector  Pour 
Location  File 
(points,  lines, 
polygons) 


MAPS  RASTERIZE 

or 

MOSS  POLYCELL 


Manual  delineation  of  watersheds  requires  preparation  of  all  the 
data  sets  from  the  conditioning  stage  and  generation  of  a  pour 
location  map.  This  map  provides  the  "seed"  or  starter  map  of 
cells  from  which  watersheds  are  delineated  by  identifying  all 
the  cells  that  flow  into  them.  The  seed  cells  represent  pour 
locations  for  watershed  pour  points  (outlets),  linear  features 
such  as  dams,  or  areas  (polygonal  features)  such  as  reservoirs 
or  multi-celled  depressions  in  the  elevation  model,  if  we  want 
to  identify  watersheds  for  the  depressions. 

This  section  discusses  the  complete  path  starting  from  genera- 
tion of  seed  points,  lines,  or  polygons  in  the  MOSS  GIS  environ- 
ment, to  development  of  watershed  maps  in  MAPS  GIS  format 
for  these  locations.  Utilize  these  processes  as  needed,  not  all 
are  necessary  for  all  projects.  This  example  illustrates  the  use  of 
the  watershed  pour  points,  however,  the  process  applies  with 
only  minor  changes  to  linear  and  polygonal  features  as  well. 


CONVERSION  PROCESS 


i 


Direct  Access 
File 


RSNAP 


Registered 
Seed  File 


WATERSHED 


Watersheds 

Generated  for 

Automatically 

Identified 
Pour  Points 
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STEP  1      Generating  a  flow  direction  map. 

The  flow  direction  data  set  developed  as  part  of  the  condition- 
ing stage  described  in  Section  II  comprises  an  essential  compo- 
nent of  manual  watershed  delineation.  Follow  the  procedures 
described  in  Part  2  of  the  HTAS  User's  Manual  to  obtain  the 
flow  direction  data  set. 


STEP  2      Generating  a  vector  pour  location  map. 

The  second  step  in  the  manual  delineation  of  watersheds 
involves  the  generation  of  a  vector  pour  location  map  in  MOSS 
using  the  GENERATE  command.  It  is  useful  to  display  flow 
accumulation  drainages  prior  to  generating  pour  locations  to 
help  place  them  on  the  drainages. 

NOTE: 

This  mafiual  assumes  advanced  knowledge  of  the  MOSS/MAPS 
CIS.  This  document  does  not  provide  examples  of  the  processing 
using  MOSS/MAPS  commands,  except  in  a  few  instances  to 
demonstrate  the  major  steps  of  the  HTAS  processing.  The  follow- 
ing example  uses  the  GENERATE  command  to  create  a  point 
map. 
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^,       EMTER  COMIIflND  ? 
I       OENERRTE 


Uhat  do  you  wish  to  call   the  new  nap  ? 


:   SEEDS. HRPU 

Enter  mapname  for  the  template  header 

<  CR=quit  ) 
:   TEHPLRTE 

ENTER  SOURCE  OF  MRP     (SOURCE  ] 

iNRHE  THE  SOURCE 

ENTER    CREATION  DATE     111/30/921 
:<CR>  OR  RNV  OTHER  DRTE  VOU  URNT  TO  PUT   IN 

ENTER    STUDVRREfl  NAME   lUOLSRG     1 
:<CR>  OR  RNV  OTHER  DESCRIPTION  VOU  URNT  TO  PUT   IN 

ENTER    DESCRIPTION 

[THE  SEED  FILE,   SRGEHL  QUAD) 

<CR>   OR  RNV  OTHER  DESCRIPTION  VOU  URNT  TO  PUT   IN 

IS  THIS  HERDER  INFORHATION  CORRECT  [V] 
:   V 

Houi  many  subjects  for  neui  map?  [      201 
:    100 

Enter  data  type  this  map  is  to  be: 

1  -  Point      2  -  Line      3  -  Polygon 
:    1 

Keyboard  or  Cursor  input  <K/C)  [CI 
:   C 

To  choose  your  point  and  enter  the  subject.     Vou  may:     Quit  <Q) 

(CR=Continue) 
:<CR> 

Enter  subject  for  this  item  (    1 
:    1 
<CURSOR   INPUT) 

To  choose  your  point  and  enter  the  subject.     Vou  may:     Quit  (Q) 

(CR=Continue) 
:<CR> 

Enter  subject  for  this  item  111 
:  2 
<CURSOR   INPUT> 

To  choose  your  point  and  enter  the  subject.     Vou  may:     Quit  <Q) 

(CR=Continue> 
:<CR> 

Enter  subject  for  this  item  [21 
:   3 
<CURSOR    INPUT> 

To  choose  your  point  and  enter  the  subject.     Vou  may:     Quit  (Q) 
<CR=Continue) 

:   Q 
****  STOP 
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STEP  3      Preparing  a  cell  (raster)  seed  map  from  the  vector  seed 
map. 

The  next  step  in  the  process  involves  preparing  a  cell  (raster) 
seed  map  from  the  vector  seed  map.  First,  RASTERIZE  the 
vector  map  into  a  type  8  raster  map,  and  then  RENUMBER  the 
values  of  all  the  zero  background  cells  to  -1. 


ENTER  COMMflMD  ? 
MRPS 

71  MAPS   IN      WORK  PROJECT  POLVGON 
130  MAPS   IN  MASTER  PROJECT  U0LSflG32 
? 

Ul    DEn.HRP  (uiindotuing  on  the  right  basemap  is  essential) 
OK    UiNDOUEO 
? 

RflSTERIZE  SEEDS. riRPU  HEIGHT  30  UIDTH  30  TVPE  8  FOR 

SEEDS. TEHPR 
PROCESSING     POINT     NUMBER  1 

PROCESSING    POINT     NUMBER  2 

PROCESSING    POINT    NUMBER  3 


PROCESSING 

POINT 

NUMBER 

79 

PROCESSING 

POINT 

NUMBER 

80 

PROCESSING 

POINT 

NUMBER 

81 

****  STOP 

OK  RASTER  1  ZED  FOR 

SEEDS. TEMPF 

RENUMBER  SEEDS. TEMPR  RS  -1 

TO  0 

FOR 

SEEDS 

OK  RENUMBERED  FOR  SEEDS. MAPR 

.MRPR 


STEP  4      Converting  MAPS  format  seed  map  into  direct  access 
format  used  in  tlie  HTAS. 

For  this  step  use  the  conversion  procedures  outlined  in  Part  1 
of  the  HTAS  User's  Manual.  For  the  purposes  of  this  descrip- 
tion, the  name  of  the  HTAS  file  into  which  the  "SEEDS.MAPR" 
has  been  converted  is  "SEEDS.MAP". 
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STEP  5      Using  RSNAP  to  correct  location  of  seed  ceiis. 


WARNING: 

The  use  of  RSNAP  applies  to  pour  point  seed  data.  It  has  not  been 
tested  with  line  or  polygon  based  seed  data.  With  these  types  of 
data  do  not  use  RSNAP  unless  you  have  tested  its  effects. 


RSNAP  finds  a  cell  with  the  nearest  flow  accumulation  value 
above  a  user-specified  flow  accumulation  threshold,  and  within 
a  threshold  distance  from  a  manually  generated  seed  cell.  The 
threshold  distance  specified  by  the  user  identifies  a  value  for 
the  number  of  cells  in  the  line  (row)  and  sample  (column) 
direction  that  represent  the  neighborhood  around  the  seed  cell. 
The  program  then  snaps  a  manually  delineated  seed  cell  to  the 
location  of  the  cell  it  found. 

The  function  of  this  program  snaps  the  manually  defined 
watershed  pour  points  or  seeds  to  the  drainage  patterns  mod- 
eled from  the  DEM.  When  manually  defining  watershed  pour 
points  there  is  the  potential  for  newly  created  seeds  not  to  be 
registered  to  the  modeled  drainage  patterns.  For  example,  a 
pour  point  generated  only  on  the  basis  of  DLG  stream  data 
might  fall  slightly  outside  of  drainage  pattern  derived  using  the 
HTAS.  In  such  a  case  the  flow  accumulation  value  for  the  point 
is  low,  because  it  is  considered  to  lie  outside  of  the  drainage 
network.  Since  few  cells  are  modeled  as  flowing  into  the  pour 
point,  the  program  delineates  only  a  very  small  watershed  for 
this  pour  point.  RSNAP  snaps  the  user  defined  pour  points  to  a 
cell  on  the  HTAS  modeled  drainage  pattern,  so  that  they  are 
properly  registered  to  the  DEM  and  full  watersheds  can  be 
delineated. 


NOTE: 

RSNAP  reads  from  and  writes  to  the  same  file,  therefore  keep  the 
original  file  in  the  HTAS  format  and  use  a  copy  for  RSNAP 
processing.  RSNAP  cannot  process  more  than  999  points. 
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OK,   R  RSNflP 
ENTER  NL,NS,FLOW  THRESHOLD, DISTANCE  THRESHOLD, SEED  FILE, COUNT 
FILE, SVSTEM(UNIX,UMS, PRIME) 

472,36Q,50, 10, 'SEEDS. HAP', 'COUNT. MAP', 'PRIME' 
nOUED  SEED  FOR    78  FROM      3,  308  TO  2,  306(2.23607  CELLS) 

MOUED  SEED  FOR    66  FROtI      9,  5110        10,  51(1.00000  CELLS) 

MOUED  SEED  FOR      2  FROM      9,  1 17  TO  8,  116(1.41421  CELLS) 

tlOUED  SEED  FOR    1 1  FROM     14,  356  TO        13,  356(1.00000  CELLS) 


MOUED  SEED  FOR    56  FROM  416,  12  TO 

MOUED  SEED  FOR    63  FROM  448,  244  TO 

MOUED  SEED  FOR    64  FROM  459,  228  TO 

TOTAL  OF  34  POINTS  MOUED 

81  SEED  CELLS  UERE  FOUND 

4c*ii<*  STOP 


417, 

12(1.00000 

CELLS) 

447, 

244(1.00000 

CELLS) 

460, 

228(1.00000 

CELLS) 

STEP  6      Delineating  waterslieds  for  tlie  user  defined  pour  locations. 

To  delineate  watersheds  for  the  generated  "seed"  map,  use  the 
"seed"  map  with  the  WTRSHED  command.  WTRSHED  finds 
watershecds  for  "seed"  cells  that  represent  user-specified  or 
automatically  assigned  pour  locations. 

To  delineate  watersheds,  "WTRSHED  utilizes  both  a  flow  direction 
data  set  and  a  "seed"  data  set.  The  "seed"  or  starter  data  set 
consists  of  background  values  of  -1  in  which  the  "seed"  cell  or 
group  of  cells  have  been  inserted  at  the  outflow  points  of  the 
desired  watersheds.  Each  start  cell  or  group  of  cells  contains  its 
own  unique  positive  value.   Using  the  flow  direction  data  set, 
WTRSHED  iteratively  reassigns  the  background  values  of  cells  in 
the  starter  data  set  to  the  value  of  the  "seed"  cell  or  group  of 
cells  into  which  the  background  cells  flow.  The  flow  direction 
map  identifies  the  flow  connectedness  and  direction  of  the 
background  cells. 


NOTE: 

Since  WTRSHED  overwrites  the  seed  file,  the  example  uses 
"SHEDS MAP",  a  copy  of  "SEEDS. MAP" for  easier  reference.  To 
keep  better  track  of  the  intermediate  products  keep  a  copy  of  the 
unchanged  seed  file,  particularly  if  disk  space  is  not  an  immedi- 
ate problem.  However,  since  you  already  have  a  version  of  it  in 
MAPS  format  and  IDIMS format,  this  is  not  necessary. 
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OK,   R  UTRSHED 

ENTER  NL,NS,DIR  FILE, MASK  FILE,SVSTEf1<UNIX,UMS,PRmE) 
472 , 369 ,  'SGML . D I R ' , 'SGML . SHED '  / PR  I  ME ' 


PASS 

PRSS 

PRSS 

PASS 

PRSS 

PASS 

PASS 

PASS 

PASS 

PASS  10 

PASS  11 

PRSS  12 

PRSS  13 

PRSS  14 


STEP  7      Converting  the  watershed  map  into  MAPS  format. 

Use  the  procedures  described  in  section  I  of  the  HTAS  Users' 
Manual.  Once  the  map  has  been  imported  into  the  MAPS  GIS,  it 
becomes  available  for  use  as  a  data  layer  in  analysis  and  carto- 
graphic modeling  using  MAPS  capabilities.  The  last  chapter  of 
this  manual  briefly  discusses  potential  applications. 
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Watersheds: 

Automatic  Delineation 


Automatic  delineation  of  watersheds  differs  from  manual  delin- 
eation in  that  the  watershed  pour  point  or  "seed"  map  is  cre- 
ated through  automatic  rather  than  manual  means.  Additionally, 
it  is  not  possible  to  generate  multi-cell  pour  locations  represent- 
ing linear  or  polygonal  features  using  the  automated  method. 


Flow 
Direction 
Data  Set 


Flow 
Accumulation 
Data  Set 


H 


DELTA 


V 


Delta 
Data  Set 


*S\-     i!?W»«»-*^' 


User  Defined 
Threshhold 


M 


SEED 


Seed 
File 


H 


WATERSHED 


Watersheds 

Generated  for 

Automatically 

Identified 
Pour  Points 
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STEP  1      Generating  the  flow  direction  and  flew  accumulatien  data 
sets  in  tlie  cenditiening  stage. 

Two  of  the  data  sets  generated  in  the  conditioning  stage  (Part 
2),  namely  the  flow  direction  data  set  and  the  flow  accumula- 
tion data  set  are  essential  for  this  procedure.  Once  these  sets 
are  generated,  create  a  delta  data  set  using  the  DELTA  com- 
mand. 


STEP  2      Identifying  flow  accumulation  increments  using  the  DELTA 
program. 

The  DELTA  program  computes  flow  accumulation  increments 
from  flow  accumulation  (COUNT)  values.  The  DELTA  process 
assists  in  automatic  delineation  of  watersheds,  by  providing  one 
of  the  data  sets  required  to  automatically  locate  watershed  pour 
points. 

Delta  value  is  the  amount  of  increase  in  flow  accumulation 
value  in  the  direction  of  flow.  Each  cell  contains  a  computed 
delta  value  derived  by  subtracting  the  flow  accumulation  value 
of  a  cell  from  the  flow  accumulation  value  of  a  cell  into  which 
it  flows.  Imagine  that  two  cells  with  a  flow  accumulation  value 
of  1000  flow  into  a  single  cell  with  a  flow  accumulation  value 
of  2000.  The  DELTA  value  of  each  of  these  cells  is  1000.  Since 
both  cells  flow  into  the  same  cell,  they  form  a  junction  of 
potential  channels.  In  the  area  where  there  is  no  juncture  of 
tributaries,  along  the  stream  channel,  the  DELTA  value  or  the 
flow  accumulation  increment  is  always  1,  because  flow  pro- 
ceeds from  one  cell  to  another  cell.  Tributary  junctures  gener- 
ally cause  a  sudden  increase  in  flow  accumulation  values. 
Therefore,  DELTA  program  generates  a  map  that  can  help  to 
identify  flow  juncture  locations,  and  the  map  created  by  it  is 
one  of  the  data  components  used  to  find  locations  of  watershed 
pour  points  using  the  SEED  program  for  watershed  delineation 
with  the  WTRSHED  program. 

The  following  example  illustrates  the  use  of  the  DELTA  com- 
mand. 

OK,  R  DELTR 

ENTER  NL,NS,DIR  FILE, COUNT  FILE, DELTA 

FILE, SVSTEM(UNIX,UnS, PRIME) 
472,369,  'DIRECT.tlflP',  'COUNT. MRP',  'DELTR. HRP',  'PRIME' 

****    3JQP 
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STEP  3      Developing  a  seed  cell  data  set  using  tlie  SEED  ppogram. 

The  next  step  of  automatic  watershed  delineation  requires 
finding  watershed  pour  points  or  seeds.  Use  the  SEED  program 
to  automatically  accomplish  this  task.  The  SEED  program  finds 
the  starter  or  "seed"  points  for  automatic  watershed  generation 
at  the  outlets  of  potential  stream  tributaries  on  the  drainage 
pattern  modeled  using  the  HTAS.  SEED  accepts  user  input  of 
the  values  (in  number  of  cells)  for  the  flow  accumulation 
threshold  above  which  the  pour  points  of  watersheds  are 
found.  It  also  uses  flow  accumulation  (COUNT)  values  and 
DELTA  values. 

For  each  cell  where  both  the  flow  accumulation  value  and  the 
delta  value  are  greater  than  the  flow  accumulation  threshold, 
SEED  assigns  to  this  cell  a  unique  positive  value.  This  means 
that,  for  a  given  threshold,  the  program  places  a  seed  where  the 
delta  value,  or  a  flow  accumulation  increment  between  two 
cells,  is  above  that  threshold,  and  where  the  flow  accumulation 
(COUNT)  value  is  above  that  threshold.  The  first  condition 
effectively  places  the  "seed"  where  there  is  a  juncture  or  tribu- 
tary relationship  between  two  potential  streams  and  the  amount 
of  flow  input  (the  increment)  from  the  streams  is  above  the 
threshold.  The  second  condition  eliminates  any  streams  which 
have  a  flow  accumulation  value  below  the  threshold  as  being 
too  insignificant  to  have  "seeds,"  thereby  eliminating  minor 
tributaries.  The  cell  that  fits  these  criteria  represents  a  pour 
point  or  a  "seed"  of  an  automatically  delineated  watershed  and 
receives  a  positive  value.  The  watersheds  grown  from  each  of 
these  pour  points  receives  the  value  of  the  "seed"  when  using 
the  WTRSHED  program.  The  background  cells  in  the  SEED  data 
set  receive  value  of  -1. 

The  following  example  illustrates  the  use  of  the  SEED  com- 
mand. 


OK,  R  SEED 

ENTER  NL,NSJHRESHOLD, EDGES, SEED  FILE, DELTA  FILE,   COUMT 

FILE, SVSTEM(UNIX,UMS, PRIME) 

472,369,3000, 1, 'SEEDS. HflP', 'DELTA. MAP', 'COUNT. HflP', 'PRIME' 

20  BRSINS  UERE  SEEDED 
****  STOP 
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STEP  4      Delineating  watersheds  fer 

tlie  automaticaiiy  derived  pour  peint  locations. 

To  delineate  watersheds  for  the  "seed"  map  generated  using  the 
SEED  program,  use  the  "seed"  map  as  one  of  the  inputs  into  the 
WTRSHED  command.  WTRSHED  finds  watersheds  for  "seed" 
cells  that  represent  automatically  assigned  pour  locations. 

To  delineate  watersheds,  WTRSHED  utilizes  both  a  flow  direc- 
tion data  set  and  the  "seed"  data  set.  The  "seed"  or  starter  data 
set  consists  of  background  values  of  -1  in  which  "seed"  cell  or 
groups  of  cells  have  been  inserted  at  the  outflow  points  of  the 
desired  watersheds,  and  each  start  cell  or  group  of  cells  has  its 
own  unique  positive  value.  Using  the  flow  direction  data  set, 
"WTRSHED  iteratively  reassigns  the  background  values  of  cells 
in  the  starter  data  set  to  the  value  of  the  "seed"  cell  or  group  of 
cells  into  which  the  background  cells  flow.  The  flow  direction 
map  identifies  the  flow  connectedness  and  direction  of  the 
background  cells. 


NOTE: 

Since  WTRSHED  overwrites  the  seed  file,  the  example  below  uses 
the  "SHEDSMAP"  name  for  easier  reference  instead  of  the  HTAS 
seed  file  "SEEDS. MAP".  To  better  track  the  intermediate  products, 
keep  a  copy  of  the  unchanged  seed  file,  particularly  if  disk  space 
is  not  an  immediate  problem.  However,  since  you  already  have  a 
version  of  it  in  MAPS  format  and  IDIMSfonnat,  this  is  not  neces- 
sary. 


OK,  R 

UTRSHED 

ENTER  NL, 

NS,DIR  FILE,riASK  FILE, 

,SYSTEri(UNIX,UMS, 

PRIME) 

472,368, 

'SGHL, 

,DIR', 

'SGHL, 

.SHED', 

'PRIME' 

PASS 

1 

PASS 

2 

PASS 

3 

PASS 

4 

PASS 

5 

PASS 

6 

PASS 

7 

PASS 

8 

PASS 

9 

PASS 

10 

PASS 

11 

PASS 

12 

PASS 

13 

' 

PASS 

14 

****  STOP 
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STEP  5      Converting  the  watershed  map  into  IVIAPS  format. 

Use  procedures  described  in  Part  1  of  the  HTAS  User's  Manual 
to  convert  watershed  maps  to  MAPS  format.  Once  the  map  is 
imported  into  the  MAPS  GIS,  it  becomes  available  as  a  data 
layer  for  analysis  and  cartographic  modeling  performed  using 
MAPS  capabilities.  The  last  chapter  of  this  manual  briefly 
discusses  these  potential  applications. 
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Delineating 

Overland  Paths  of  How 


To  delineate  overland  flow  paths  for  various  point  sources  on 
the  landscape  use  a  procedure  similar  to  the  procedure  for 
manual  delineation  of  watersheds.  First,  generate  a  starter  or 
"seed"  data  set  which  in  this  case  does  not  represent  watershed 
pour  points,  but  sources  of  flow.  These  sources  of  flow  include 
springs  or  point  sources  of  pollution.  The  HTAS  OVERLAND 
program  delineates  the  steepest  path  of  flow  along  the  land- 
scape using  the  flow  direction  data  set,  by  tracing  out  the  cells 
into  which  the  initial  cell  flows. 


STEP  1      Generating  a  flow  direction  map. 

The  flow  direction  data  set  developed  as  part  of  the  condition- 
ing stage  described  in  Part  2  is  an  essential  component  of 
automatic  watershed  delineation.  Follow  the  procedures  de- 
scribed in  Part  2  of  the  HTAS  User's  Manual. 


STEP  2      Generating  a  vector  point  source  map. 

The  second  step  of  the  automatic  delineation  of  watersheds 
involves  generating  a  vector  point  source  map  in  MOSS  using 
the  GENERATE  command.  This  command  allows  the  user  to 
generate  the  points  anywhere  on  the  landscape. 

NOTE: 

This  manual  assumes  advanced  knowledge  of  the  MOSS/MAPS 
GIS,  and  examples  of  the  processing  using  MOSS/MAPS  com- 
mands will  not  be  generally  provided.  However,  a  few  exceptions 
will  be  made  so  as  to  demonstrate  the  major  steps  of  the  HTAS 
processing.  An  example  of  using  the  GENERATE  command  to 
create  a  point  map  is  as  follows: 
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ENTER  COMMflHD  ? 
GENERRTE 
Uhat  do  you  wish  to  call   the  neui  map  ? 


:   SOURCES. riRPU 
Enter  mapname  for  the  template  header 

(  CR=quit   ) 

TEHPLflTE 
ENTER  SOURCE  OF  NflP     [SOURCE  ] 

NRHE  THE  SOURCE 
ENTER    CREATION  DATE     [11/30/921 
<CR>  OR  RMV  OTHER  DRTE  VOU  URNT  TO  PUT   IN 
ENTER    STUOVAREA  NAME   [UOLSAG     1 

<CR>   OR  RNV  OTHER  DESCRIPTION  VOU  URNT  TO  PUT    IN 
ENTER    DESCRIPTION 
[THE  POINT  SOURCES  FILE,   SAGEHL  QUAD! 

<CR>  OR  RNV  OTHER  DESCRIPTION  VOU  URNT  TO  PUT  IN 

IS  THIS  HERDER  INFORNATION  CORRECT  [Vl 
:  V 

Hou)  many  subjects  for  new  map?  [   201 
:  100 

Enter  data  type  this  map  is  to  be: 

1  -  Point   2  -  Line   3  -  Polygon 
:  1 

Keyboard  or  Cursor  input  (K/C>  [CI 
:  C 

To  choose  your  point  and  enter  the  subject 

(CR=Continue) 
:<CR> 

Enter  subject  for  this  item  [    ] 
:    1 
<CURSOR    INPUT) 

To  choose  your  point  and  enter  the  subject 

<CR=Continue) 
:<CR> 

Enter  subject  for  this  Item  [1    ] 


Vou  may:     Quit  (Q) 


Vou  may:     Quit  <Q) 


<CURSOR    INPUT> 

To  choose  your  point  and  enter  the  subject 

(CR=Continue) 
:<CR> 

Enter  subject  for  this  item  [21 
:   3 
<CURSOR    INPUT> 

To  choose  your  point  and  enter  the  subject 

(CR=Continue> 
:   Q 
****  STOP 


Vou  may:     Quit  (Q) 


Vou  may:     Quit  <Q) 
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STEP  3      Preparing  a  cell  (raster)  seed  map  from  the  vector  seed 
map. 

The  next  step  prepares  a  cell  (raster)  point  source  seed  map 
from  the  vector  point  source  seed  map.  First  RASTERIZE  the 
vector  map  into  a  type  8  raster  map,  and  then  RENUMBER  the 
values  of  all  the  zero  background  cells  to  -1. 


ENTER  COMMAND  ? 
MRPS 

71  MAPS   IN      WORK  PROJECT  POLVGON 
130  MAPS   IN  MASTER  PROJECT  U0LSAG32 
? 

Ul    DEn.MRP  (uiindoiuing  on  the  right  basemap  is  essential) 
OK     UINDOUED 
? 

RASTERIZE  SOURCES. riRPU  HE  30  Ul    30  TVPE  6  FOR 
SOURCES. TEHPR 

PROCESSING    POINT    NUMBER  1 

PROCESSING    POINT    NUMBER  2 

PROCESSING    POINT    NUMBER  3 


PROCESSING 

POINT 

NUMBER 

79 

PROCESSING 

POINT 

NUMBER 

80 

PROCESSING 

POINT 

NUMBER 

81 

****  STOP 

OK  RASTER  1  ZED  FOR 

SOURCES 

.TEMPR 

RENUMBER  SOURCES, 

,  TEHPR  ( 

IS  -1 

TO 

OK  RENUMBERED  FOR  SOURCES 

.MRPA 

0  FOR  SOURCES. HRPR 


STEP  4      Converting  MAPS  format  point  source  map  into  direct 
access  format  used  in  tiie  HTAS. 

For  this  step,  use  the  conversion  procedures  outlined  in  Part  1 
of  the  HTAS  User's  Manual.  For  the  purposes  of  this  descrip- 
tion, the  HTAS  file,  SOURCES.MAPR  was  converted  to  an  HTAS 
file  named  SOURCES.MAP. 
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STEP  5      Delineating  overland  paths  for  the  user  defined  point 
sources. 

To  delineate  overland  for  the  generated  point  source  map,  use 
the  point  source  map  with  the  OVERLAND  program.   OVER- 
LAND traces  overland  paths  for  user-specified  cell(s)  contained 
in  a  starter  point  sources  data  set.  The  program  uses  the  flow 
direction  data  set  to  produce  a  path  or  paths  by  following  the 
cell-to-cell  linkage  (how  cells  flow  into  each  other)  starting  from 
the  user  specified  cell(s)  until  it  reaches  the  data  set  edge.    The 
results  of  this  process  are  helpful  in  tracking  the  overland  path 
from  point  source  of  pollution  into  nearby  streams. 


OK,  R  OUERLRND 

ENTER  NL,NS,DIR  FILE, MASK  FILE,SVSTEn(UNIX,UMS,PRmE) 
472,368,  'DIRECT. nflP',  'SOURCES. HRP', 'PRIME' 

PASS  1 

PASS  2 

PASS  3 

PASS  4 

:tf***    STOP 


STEP  6      Converting  the  overland  paths  map  into  MAPS  format. 

Again,  use  procedures  described  in  Part  1  of  the  HTAS  User's 
Manual.  Once  the  map  is  imported  into  the  MAPS  GIS,  it  be- 
comes available  for  use  as  a  data  layer  in  analysis  and  carto- 
graphic modeling  performed  using  MAPS  capabilities.  The  last 
chapter  of  this  manual  briefly  discusses  potential  applications. 
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Locating  Bevation  Minima  & 

Delineating  Drainage  Patterns 


There  two  ways  to  delineate  drainage  patterns  using  the  HTAS. 
The  first  method  uses  the  flow  accumulation  data  set.  By 
displaying  the  values  of  the  flow  accumulation  map  developed 
in  the  Conditioning  Stage  at  various  increments,  for  example 
50-100,  100-999,  999-50000,  a  drainage  pattern  appears  on  the 
screen.  The  display  of  the  drainage  network  gradually  increases 
in  complexity  as  the  bottom  threshold  of  flow  accumulation 
values  being  displayed  decrease.  Using  MAPS  commands  such 
as  EXTRACT,  extract  ranges  of  values  to  obtain  a  network  with 
a  density  that  approximates  the  level  of  detail  desired,  for 
example  approximating  the  density  of  a  stream  network  derived 
from  DLG  data. 

Use  the  programs  MINIM  and  OVERLAND  together  to  obtain  a 
map  of  drainage  patterns  via  an  alternative  method.  MINIM 
finds  all  the  local  elevation  minima,  or  downward  inflection 
points  in  a  depressionless  DEM.  OVERLAND  uses  the  data  set 
produced  by  MINIM  to  delineate  drainage  patterns  by  connect- 
ing the  local  minima  with  flow  paths.  According  to  Jenson,  this 
method,  as  opposed  to  the  flow  accumulation  method,  pro- 
duces narrow  drainage  paths  for  areas  of  broad,  shallow  valleys 
where  the  flow  accumulation  data  set  would  produce  a  drain- 
age network  with  wide  lines. 
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STEP  1      Locating  elevation  minima. 

To  locate  elevation  minima,  such  as  downward  inflection  points 
and  flats,  use  the  HTAS  program  MINIM.  The  MINIM  program 
reads  an  elevation  file  and  examines  each  cells  neighborhood 
within  a  3x3  window  (a  nine  cell  neighborhood).  If  any  of  the 
endpoints  of  symmetric  cross-sections  through  the  neighbor- 
hood are  higher  than  the  center,  then  the  center  cell  receives  a 
value  of  2  to  represent  the  bottom  of  a  valley,  or  a  V-shape  on 
the  landscape.  If  the  cell  is  equal  in  value  to  all  its  neighbors,  it 
is  a  flat  and  receives  the  value  of  3.  MINIM  delineated  cells,  by 
identifying  valley  shaped  areas  and  flats,  form  a  crude 
unconnected  drainage  network. 


OK,  R  nmin 

ENTER  NL^NSJOPO  FILE^MmiMfl  FILE,SVSTEM(UNIX,UMS,PRmE) 
472, 360 ,  'DEM . nflP ' , 'n I M I nfl . HflP ' , ' PR  I  HE • 
COMPLETED  fl  MOREFLflTS  PASS 


STEP  2      Delineating  a  drainage  network. 

Running  the  HTAS  program  OVERLAND  on  the  MINIMA  data 
set  generates  a  map  in  which  the  overland  paths  for  all  the 
MINIMA  cells  are  traced  using  the  flow  direction  data  set.  This 
produces  an  fully  connected  drainage  network  data  set. 

The  OVERLAND  program  traces  overland  paths  for  a  starter 
data  set  consisting  of  cells  assigned  positive  values  in  a  back- 
ground of  -1.  In  this  case,  the  starter  data  set  is  the  MINIM  data 
set.  The  program  uses  the  flow  direction  map  to  produce  paths 
by  following  the  cell-to-cell  linkage  (how  cells  flow  into  each 
other)  starting  from  the  user  specified  cell(s)  until  the  program 
reaches  the  data  set  edge. 


OK,   R  OUERLRND 

ENTER  NL,NS,DIR  FILE, MASK  FILE,SVSTEN<UNIX,UMS,PRmE) 
472, 36fl , 'D I RECT . HAP ' , 'H I N I H . HflP ' , 'PR  I  HE ' 

PASS  1 

PASS  2 

PASS  3 

PASS  4 
****  STOP 
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Generating  a  Table  of 

lAbtershed  Pour  Points 


Create  a  table  of  pour  points  using  the  HTAS  program  PPT  ABLE 
to  review  the  location  of  watershed  pour  points  or  find  out 
their  elevation.  PPTABLE  builds  a  table  of  pour  points  for 
watersheds.  To  accommodate  data  sets  where  there  are  internal 
or  closed  drainage  patterns  and  where  pour  points  are  very 
close  or  equal  in  elevation  for  a  watershed  A  to  drain  to  either 
watershed  B  or  C,  the  program  computes  all  possible  pour 
points.  The  table  provides  information  on  the  label,  elevation 
and  location  of  the  pour  point  on  the  map  in  line  (row)  and 
sample  (column)  format. 

OK,  R  PPTflBLE 
ENTER  NL,NS,BflSIM  FILE,ELEU  FILE, TABLE 

FILE, SVSTEn(UNIX,Uf1S, PRIME) 
472, 36Q, 'SEED. MRP', 'DEM. MRP', 'PPTRBLE. SEED', 'PRIME' 

Example  of  output: 

OK,  SLIST  PPTRBLE.SEED 

1  THIS   IS  THE  COMPLETE  LIST  OF  POUR  POINTS 


SEQUENCE  NO. 

LABEL 

PAIR 

ELEURTION 

LINE 

SAMP 

1 

0/ 

1 

2070 

2 

0 

2 

0/ 

2 

2151 

4 

0 

3 

0/ 

3 

2146 

7 

0 

4 

0/ 

4 

2126 

7 

0 

163 

55/ 

56 

2608 

448 

49 

164 

0/ 

50 

2485 

453 

0 

165 

0/ 

56 

2381 

469 

0 

166 

0/ 

54 

2194 

463 

0 
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HTAS  Commands 

Reference 
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This  part  describes  each  of  the  commands  used  by  HTAS, 
listing  them  in  alphabetical  order. 

CONCAVE  -  Fills  large  depressions  in  a  DEM.  CONCAVE 
requires  specifications  for  the  size  of  box  or  window 
(BOXNL,  BOXNS)  and  the  dimensions  by  which  the  boxes 
overlap  in  both  line  (row)  and  sample  (column)  by  a  value 
specified  as  OVERLAP.  The  largest  diameter  of  the  largest 
depression  to  be  filled  must  be  less  than  the  overlap.  The  size 
of  the  box  (BOXNL,  BOXNS)  must  not  be  greater  than  400. 
CONCAVE  is  a  part  of  the  conditioning  stage  for  delineating 
watersheds,  however  its  product  also  has  value  in  identifying 
surface  depressions  in  the  DEM.  CONCAVE  reads  from  and 
writes  to  the  original  DEM  file,  so  it's  good  to  keep  a  copy  of 
the  original.  The  following  example  uses  a  copy  of  the  DEM 
file  called  CONCAVE.MAP. 

Input  parameters: 

Number  of  lines  or  rows  (NL), 

Number  of  samples  or  columns  (NS), 

Window  size  in  lines  or  rows  (BOXNL), 

Window  size  in  samples  or  columns  (BOXNS), 

Overlap  between  windows  (OVERLAP), 

Name  of  the  elevation  file  in  single  quotes  (LEV  FILE), 

Name  of  the  system  in  single  quotes,  e.g.,  'PRIME'  (SYSTEM) 

OK,   R  CONCRUE 

ENTER  NL,NS, BOXNL, BOXMS,OUERLflP, LEU  FILE,SVSTEt1<UNIX,UMS,PRmE) 
472, 360 ,' 100, 100, 60, ' CONCRUE . MRP ' , ' PR  I  ME ' 
System  =  PRIME 

Open  as  a  Prime  direct  access  file,  ns  =     369 
SGHL.CONCflUE 
opened  successful ly 


SL  = 

1 

SL  = 

41 

SL  = 

81 

SL  = 

121 

SL  = 

161 

SL  = 

201 

SL  = 

241 

SL  = 

281 

SL  = 

321 

SL  = 

361 

SL  = 

401 

SL  = 

441 

****  STOP 
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CONVERT  -  Converts  from  and  to  sequential  access  raster  transfer 
map  file  (IDIMS  format)  and  a  direct  access  raster  map  file  used 
by  the  HTAS. 

Input  parameters: 

Select  transfer-to-direct  (option  0)  or  direct-to-transfer  (option  1) 

If  direct-to-transfer  select  1*2  datatype  (option  2)  or  1*4  datatype 
(option  4) 

Number  of  lines  (rows) 

Number  of  samples  (columns) 

Converting  from  sequential  access  raster  transfer  format  to  direct 
access  format: 

OK,  R  CONUERT 

Conversion  desired  (0=1*2  transfer-to-direct,    1=1*2  or 
1*4  direct-to-transfer): 
0 

Transfer  DEtl  input  image  name: 

DEn.nRP.i6 

Direct  access  output  DEM  image: 
DEn.HflP 

Lines  in   image: 
472 

Samples  in  image: 

Converted  200  I ines 

Converted  400  lines 

mil**  STOP 

Converting  from  direct  to  sequential  access  raster  transfer  format: 

OK,  R  CONUERT 

Conversion  desired  (0=1*2  transfer-to-direct,  1=1*2  or 
1*4  direct-to-transfer): 
1 

Direct  access  input  image  name: 
SHED.nRP 

Transfer  output  image  name: 
SHED.nRP. 16 

Direct-to-transfer  datatype  (2=1*2,  4=1*4): 
2 

Lines  in  image: 
472 

Samples  in  image: 
369 

Converted        200  I ines 

Converted        400  I ines 
****   STOP 


NOTE: 

To  CONVERT  1*4  files,  select  option  4  in  "Direct-to-transfer' 
datatype. 
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COUNT  -  Computes  flow  accumulation  values  from  flow  direction 
values.  COUNT  assigns  each  cell  in  the  COUNT  data  set  a  value 
equal  to  the  number  of  cells  that  flow  into  it.  Cells  having  the  flow 
accumulation  value  of  zero  (to  which  no  other  cells  flow)  generally 
correspond  to  the  pattern  of  ridges.  Because  all  cells  in  a 
depressionless  DEM  (a  product  of  FILLSNGL  and/or  CONCAVE) 
contain  a  path  to  the  data  set  edge,  the  pattern  formed  by  highlight- 
ing cells  with  values  higher  than  some  threshold  value  delineates  a 
fully  connected  drainage  network.  As  the  threshold  value  decreases, 
the  density  of  the  drainage  network  increases.  COUNT  is  a  part  of 
the  conditioning  stage  for  delineating  watersheds,  however  its  prod- 
uct also  has  a  value  for  delineating  drainage  networks  and  potential 
stream  flow  modeling. 

Input  parameters: 

Number  of  lines  or  rows  (NL), 

Number  of  samples  or  columns  (NS), 

Restart  option  if  processing  interrupted  (CODE  1  TO  RESTART), 

Name  of  flow  direction  file  in  single  quotes  (DIR  FILE), 

Name  of  the  flow  accumulation  file  in  single  quotes  (COUNT  FILE), 

Name  of  the  system  in  single  quotes,  here  it  is  'PRIME'  (SYSTEM) 

OK,  R  COUNT 

ENTER  NL,NS,CODE   1  TO  RESTART, DIR  FILE, COUNT  FILE,SVSTEM<UNIX,UriS,PRmE) 
472,369,0,  'DIRECT. nflP', 'COUNT. MAP', 'PRIME' 

DOWNURRD  PRSS    1 

UPUflRD  PASS    2 

DOUNUARD  PASS    3 

UPURRD  PASS    4 

DOUNUARD  PRSS    5 

UPUARD  PASS    6 

DOUNUARD  PASS    7 

UPURRD  PASS    8 

DOUNUARD  PASS    9 

UPUARD  PASS   10 

DOUNUARD  PASS   11 

UPURRD  PRSS   12 

DOUNUARD  PASS   13 

UPUARD  PRSS   14 
in***  9T0P 
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DELTA  -  Computes  flow  accumulation  increments  from  flow 
accumulation  (COUNT)  values.  DELTA  is  a  process  used  in 
automatic  delineation  of  watersheds  where  the  user  provides  a 
threshold  flow  accumulation  value  at  which  watershed  pour 
points  should  be  inserted.   DELTA  value  represents  the  amount  of 
increase  in  flow  accumulation  value  in  the  direction  of  flow. 
The  program  computes  DELTA  values  for  every  cell  by  subtracting 
the  flow  accumulation  value  of  a  cell  from  the  flow  accumulation 
value  of  a  cell  into  which  it  flows.  Imagine  that  two  cells  with  a 
flow  accumulation  value  of  1000  flow  into  a  single  cell  with  a  flow 
accumulation  value  of  2000.  The  DELTA  value  of  each  of  these 
cells  is  1000.  Since  they  both  flow  into  the  same  cell,  they  form  a 
junction  of  potential  streams  or  channels.  In  the  area  where  there  is 
no  juncture  of  tributaries  along  the  stream  channel,  the  DELTA 
value  or  the  flow  accumulation  increment  is  always  1  because 
flow  proceeds  from  one  cell  to  another  cell.  At  the  juncture  of 
tributaries  a  sudden  increase  in  flow  accumulation  values  is 
realized.  Therefore,  DELTA  program  identifies  locations  of  flow 
junctures,  and  the  map  created  by  it  is  one  of  the  data  components 
used  to  find  locations  of  watershed  pour  points  using  the  SEED 
program  for  watershed  delineation  with  the  WTRSHED  program. 
(See  SEED  and  WTRSHED). 

Input  parameters: 

Number  of  lines  or  rows  (NL), 

Number  of  samples  or  columns  (NS), 

Name  of  flow  direction  file  in  single  quotes  (DIR  FILE), 

Name  of  the  flow  accumulation  file  in  single  quotes  (COUNT  FILE), 

Name  of  the  delta  file  in  single  quotes  (DELTA  FILE), 

Name  of  the  system  in  single  quotes,  here  it  is  'PRIME'  (SYSTEM) 

OK,    R  DELTR 

ENTER  NL,NS,DIR  FILE, COUNT  FILE, DELTA  FILE,SVSTEn(UNIX,UMS,PRmE) 
472, 365 /D I RECT .  HRP ' ,  'COUNT .  HflP ' ,  'DELTR .  MRP '  /PR  I  HE ' 
*,ti**  STOP 


HTAS  REFERENCE 


Program  Listings [61 


DIRECT  -  Computes  flow  direction  values  from  a  DEM.  The 
flow  direction  for  a  cell  is  the  direction  water  flows  out  of  the 
cell,  and  the  flow  direction  value  corresponds  to  the  orienta- 
tion of  one  of 
the  eight  cells  that  surround  it  as  follows: 

64(NW)  128(N)  1(NE) 

32(W)  center  2(E) 

16(SW)  8(S)  4(SE) 

There  are  four  possible  conditions  to  consider  in  determining 
flow  direction.  The  first  condition  applies  only  to  single  cell 
depressions,  when  the  immediate  surrounding  cells  have  a  value 
higher  than  the  center  cell,  a  condition  true  only  in  the  unfilled 
DEMs.  In  that  case  DIRECT  assigns  a  negative  value  to  the 
depression  and  provides  a  message  stating  "COULD  NOT  SOLVE 
FOR  ALL  CELLS".  The  second  condition  occurs  when  the  dis- 
tance-weighted drop  from  the  center  cell  is  higher  for  one  of  the 
surrounding  cells  than  for  the  other  neighborhood  cells.  A  flow 
direction  value  is  assigned  to  the  center  cell  that  leads  to  the  cell 
where  the  drop  (or  the  slope  gradient)  is  the  steepest.  Most  of 
the  situations  are  like  this.  In  the  third  case,  there  are  two  or 
more  cells  that  have  the  same  steepest  drop,  and  one  has  to  be 
selected  as  the  direction  of  flow.  This  is  accomplished  using  a 
logical  table  look-up  operation,  where,  for  example,  if  three 
adjacent  cells  have  the  same  steepest  drop  from  the  center  cell, 
the  program  selects  the  middle  one  of  these  cells  for  flow 
direction.  Another  possibility  for  this  example  is  if  two  cells  on 
opposite  sides  have  an  equal  drop,  then  the  program  arbitrarily 
choses  one  of  them.  In  the  fourth  condition  the  center  cell  is 
located  in  a  flat  area  and  all  cells  are  equal  (  or  greater)  in 
elevation.  The  outflow  point  is  not  known.  Consequently, 
processing  requires  that  all  the  cells  belonging  to  the  first, 
second,  or  third  condition  be  resolved,  and  then,  in  an  iterative 
process,  assigns  flow  direction  to  neighboring  cells  if  the  neigh- 
bor has  a  defined  flow  direction  that  does  not  point  back  to  the 
tested  cell.  In  this  way  the  flow  direction  assignments  grow  into 
the  flat  area  from  the  flats'  outflow  point  until  all  the  cells  have 
flow  directions  assigned.  DIRECT  is  a  part  of  the  conditioning 
stage  for  delineating  watersheds,  however  its  product  also  has  a 
value  for  helping  visualize  general  flow  patterns  on  the  land- 
scape. It  can  be  applied  to  an  unfilled  DEM  or  to  a  filled  DEM. 
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Input  parameters: 

Number  of  lines  or  rows  (NL), 

Number  of  samples  or  columns  (NS), 

Restart  option  if  processing  interrupted  (CODE  1  TO  RESTART), 

Name  of  the  elevation  file  in  single  quotes  (LEV  FILE), 

Name  of  the  flow  direction  file  in  single  quotes  (DIR  FILE), 

Name  of  the  system  in  single  quotes,  here  it  is  'PRIME'  (SYSTEM) 

OK,  R  DIRECT 
ENTER  NL,NS,CODE   1  TO  RESTART, LEU  FILE, DIR  FILE,SVSTEM(U 

MIX,UMS,PRmE) 
472,36g,0,  'DEn.nflP',  'DIRECT. nflP'/PRIHE' 
Open  ELEU  and  DIR  f i les,   ns  =  359 

DEM. MRP 

opened  successful ly 
DIRECT. MRP 

opened  successful ly 
DOUNURRD  PASS  1FIRSTL  2LflSTL  471 
UPURRD  PASS   2FIRSTL  7LASTL  465 
DOUNUARD  PASS  '3FIRSTL  7LflSTL  461 
COULD  NOT  SOLUE  FOR  RLL  CELLS 
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FILLSNGL  -  Fills  single-cell  depressions  in  a  DEM.  The  program 
looks  at  each  cell  in  comparison  to  its  adjacent  eight  cells  in  a 
3x3  window  (i.e.,  a  nine-cell  neighborhood).  If  the  cell  is 
lower  in  value  than  all  of  its  neighbors,  the  program  raises  the 
cell's  value  to  that  of  its  lowest  neighbor.   While  CONCAVE 
also  fills  single  cell  depressions,  this  program  generates  a 
data  set  which  helps  in  distinguishing  between  single-cell  and 
multi-cell  depressions.  FILLSNGL  reads  from  and  writes  to  the 
original  DEM  file,  so  make  a  copy  of  the  original. 

Input  parameters: 

Number  of  lines  or  rows  (NL), 

Number  of  samples  or  columns  (NS), 

Name  of  the  unprocessed,  original  elevation  file  in  single 
quotes  (TOPO  FILE), 

Name  of  the  system  in  single  quotes,  here  it  is  'PRIME'  (SYS- 
TEM) 


OK,  R  FILLSNGL 

ENTER  NL,NS,LEU  FILE,SVSTEM<UNIX,UMS,PRmE) 
472,368, 'DEM. tlflP'/PRiriE' 
Open  as  a  Prime  direct  access  file,  ns  = 

DEM. MAP 

opened  successful ly 
^ii**   5JQP 


369 
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IMPORT  (MAPS)  -  this  MAPS  command  permits  conversion  and 
import  of  sequential  access  raster  transfer  (IDIMS)  files  into 
the  MAPS  system.  The  Moss  User's  Manual  describes  the 
IMPORT  command  in  detail.  What  follows  is  a  brief  descrip- 
tion of  new  IMPORT  capabilities  which  permit  import  of  1*4 
maps  into  MAPS.  Install  the  new  IMPORT  version  included 
with  the  HTAS  package  in  a  directory  of  your  choice  and  set 
up  PTHBLK.DT  in  the  directory  from  which  you  are  running 
MOSS/MAPS  so  that  when  IMPORT  is  used  the  new  version  is 
found.  The  new  version  is  provided  in  files  IDIMS. RUN  and 
IDIMPR.RUN. 


Below  is  an  example  of  how  the  MAPS  IMPORT  command 
works  with  the  1*4  maps.  The  import  1*4  maps  use  the 
"DOUBLEWORD"  instead  of  the  "WORD"  option. 

OK,  R  IDIHS 

SINGLE  COMMAND  UERSION  OF  MAPS 
ONLV  BVE  AND  IMPORT  WILL  WORK 

0  MAPS  IN    WORK  PROJECT  POLVGON 
NO  MASTER  PROJECT  IS  OPEN 
? 

IMPORT  COUNT. MRP.  16  FORMRT  IDIMS  DOUBLEWORD  FOR  COUNT. MRP 
TVPE  8 

ENTER  THE  FOLLOWING  INFORMATION  FOR  INPUT  MAP 

NUMBER  OF  ROWS  ? 
472 

NUMBER  OF  COLUMNS  ? 
369 

CELL  HEIGHT  <IN  METERS)  ? 
30 

CELL  WIDTH  <IN  METERS)  ? 
30 

ENTER  THE  MINIMUM  BOUNDING  RECTANGLE 

MINIMUM  X  (WEST)  ? 
702750 

MAXIMUM  X  (EAST)  ? 
713820 

MINIMUM  V  (SOUTH)  ? 
4416390 

MAXIMUM  V  (NORTH)  ? 
4430550 

ENTER  MAP  PROJECTION  INFORMATION 

PROJECT  1 0N(0-20)  ? 
1 

ELL  IPSO  I D(0- 19)  ? 
0 

LONGITUDE  OF  ANV  POINT  WITHIN  THE  UTM  ZONE  ? 
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LATITUDE  OF  flMV  POINT  Nl THIN  UTM  ZONE       12  ? 
35 

ENTER  MAP  DESCRIPTION  INFORMATION 

HAP  DESCRIPTION  ? 
URTERSHEDS  BRSED  ON  THRESHOLD  OF  3000 

STUDV  AREA  ? 
UOLSRG 

PROJECTION  DESCRIPTION  ? 
UTM 

DATE  ? 
10/01/92 

SOURCE  ? 
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MAPSTOIDIMS  -  Permits  conversion  of  raster  files  in  MAPS 
format  into  the  sequential  access  raster  transfer  format 
(IDIMS)  as  an  operating  system  (PRIMOS)  level  utility.  Origi- 
nally written  to  convert  and  write  files  to  tape,  this  command 
also  converts  and  writes  files  to  disk  by  using  defaults. 
MAPSTOIDIMS  creates  files  with  .16  and  .TH  suffixes  after  the 
name  of  the  input  file.  The  .16  files  are  the  raster  files  used 
for  processing,  while  the  .TH  files  contain  information  about 
various  parameters  of  the  file,  such  as  number  of  rows  and 
columns  or  coordinates  of  the  Minimum  Bounding  Rectangle 
(the  window  of  the  map)  to  which  the  user  may  refer. 

OK,  tIRPSTOIDinS 

MAPS  TO   IDIMS  CONVERSION 

Enter  macro  file  name  to  execute  file  copy 

(  Default  is  TMCOPV   > 
<CR> 

Enter  tape  drive  name  MTO  or  MT1   (  Default  =  MTO  ) 
<CR> 

Enter  tape  drive  density  (800,    1600  or  6250;  Default  =  1600  )  : 

How  many  maps  are  already  on  this  tape  ? 
:<CR> 

Please  enter  the  DIRECTORV  (under  NOSSOflTfl)  or  WORK 

and  the  .DT  (project  master  or  POLVGON   : 

i.e.,  UOLSflG  UOLSflG.DT    or    WORK  POLVGON. DT 
:   UOLSRO  32  U0LSRG32.DT 
Enter  map  name  (  CR  =  End  ) 
:    DEnFILE 

Map-  DEMFILE      has  been  processed. 
Enter  map  name  (  CR  =  End  ) 
:<CR> 
*   INSTRUCTIONS  TO  COMPLETE  MflPSTOIDIMS  * 

1)  List  out  the  .TH  file  created  In  the 

same  directory  that  the  MflPSTOIDIMS  was  i 

executed  (SPOOL  or  SLIST  on  the  PRIME).  ' 

This  file  contains  the  MBR,  cell  size, 

and  projection  criteria  of  the  file, 

and  needs  to  accompany  the  TflPE  , 

trans fered  to  IDIMS  | 

2)  R  the  macro  called  TMCOPV  to  transfer  ( 
the  cell  file  to  magnetic  tape:  j 
•  Mount  the  tape 
t                 R  TMCOPV  (CPL  TMCOPV) 

MflPSTOIDIMS  completed. 
«)icitt>ic  STOP 

WARNING: 

MAPSTOIDIMS  works  only  with  16-bit  (word  length)  integer  data 
and  some  of  the  MAPS  real  number  files  32  (double  word)  data. 
To  correct  tbisprobletn  use  the  FUNCTION  command  INTEGER 
option  in  maps  which  will  convert  the  real  number  32-bit  map 
into  an  integer  16-bit  map. 
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MEVIM  -  Finds  local  minima  and  flats  in  a  DEM.  The  MINIM 
program  reads  an  elevation  file  and  examines  each  cells 
neighborhood  within  a  3x3  window  (a  nine  cell  neighbor- 
hood). If  any  of  the  endpoints  of  symmetric  cross-sections 
through  the  neighborhood  are  higher  than  the  center,  then 
the  program  gives  the  center  cell  a  value  of  2  to  represent  the 
bottom  of  a  valley,  or  a  V-shape  on  the  landscape.  If  the  cell 
is  equal  in  value  to  all  its  neighbors,  it  is  a  flat  and  receives  a 
value  of  3.  MINIM-delineated  cells,  which  identify  valley- 
shaped  areas  and  flats,  form  a  crude  unconnected  drainage 
network.  Use  OVERLAND  to  connect  these  cells  and  follow 
the  flow  path  between  them  to  delineate  a  drainage  pattern. 

Input  parameters: 

Number  of  lines  or  rows  (NL), 

Number  of  samples  or  columns  (NS), 

Name  of  the  unprocessed,  original  elevation  file  in  single 
quotes  (TOPO  FILE), 

Name  of  the  new  minima  file  in  single  quotes  (MINIMA  FILE) 

Name  of  the  system  in  single  quotes,  here  it  is  'PRIME' 
(SYSTEM) 

OK,  R  niNin 

ENTER  NL,NS,TOPO  FILE.MINmfl  FILE,SVSTEH(UNIX,UMS,PRmE) 
472,368,  'DEn.tlflP',  'HINinfl. MAP',  'PRIME' 

COMPLETED  fl  MOREFLflTS  PASS 
****  STOP 


Program  Listings 


eitiHTASReerence 


OVERLAND  -  Traces  overland  paths  for  user-specified  cell(s) 
contained  in  a  starter  data  set.  The  program  uses  the  flow 
direction  data  set  to  produce  a  path  or  paths  by  following  the 
cell-to-cell  linkage  (how  cells  flow  into  each  other)  starting 
from  the  user  specified  cell(s)  until  it  reaches  the  data  set 
edge.  Create  the  starter  data  set  by  generating  the  points  in 
MOSS  to  represent  point  sources,  rasterizing  them  to  MAPS 
format,  converting  them  using  MAPSTOIDIMS  and  CONVERT 
to  the  direct  access  file  format  used  by  the  HTAS.  The  OVER- 
LAND process  assigns  the  starter,  or  point  source  cell  label  to 
all  the  cells  it  flows  into,  i.e.,  to  the  entire  path.  The  results  of 
this  process  are  helpful  in  tracking  the  overland  path  from 
point  sources  of  pollution  into  nearby  streams.  When  using 
OVERLAND  on  a  data  set  created  by  MINIMA,  it  delineates  a 
complete  drainage  network  pattern  for  the  modeled  area. 

Input  parameters: 

Number  of  lines  or  rows  (NL), 

Number  of  samples  or  columns  (NS), 

Name  of  the  flow  direction  file  in  single  quotes  (DIR  FILE), 

File  containing  points  for  which  the  path  is  modeled 
(MASK  FILE), 

Name  of  the  system  in  single  quotes,  here  it  is  'PRIME' 
(SYSTEM) 

OK,  R  OUERLRND 

ENTER  NL,NS,DIR  FILE, MASK  FILE,SVSTEM<UNIX,UriS,PRmE) 
472, 36Q,  'DIRECT. HflP',  'SOURCES. MRP',  'PRIME' 

PASS  1 

PASS  2 

PASS  3 

PASS  4 

****  STOP 
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PPTABLE  -  Builds  a  table  of  pour  points  for  watersheds. 
To  accomodate  data  sets  where  there  are  internal  or  closed 
drainages  and  where  pour  points  are  very  close  or  equal  in 
elevation  for  watershed  A  to  drain  to  either  watershed  B  or  C, 
this  program  computes  all  possible  pour  points.  For  each 
watershed,  the  program  examines  every  bordering  watershed 
(including  the  zero  or  off  the  edge  watersheds).  The  table 
includes  information  on  the  label,  elevation  and  location  of 
the  pour  point  on  the  map  in  line  (row)  and  sample  (column) 
format. 

Input  parameters: 

Number  of  lines  or  rows  (NL), 

Number  of  samples  or  columns  (NS), 

File  of  watersheds  delineated  using  WTRSHED  in  single 
quotes  (BASIN  FILE) 

Name  of  the  elevation  file  in  single  quotes  (LEV  FILE), 

Name  of  the  table  to  be  produced  in  single  quotes 
(TABLE  NAME) 

Name  of  the  system  in  single  quotes,  here  it  is  'PRIME' 
(SYSTEM) 

OK,   R  PPTflBLE 
ENTER  NL,NS,BflSIN  FILE,ELEy  FILEJflBLE 

FILE, SVSTEM(UNIX,UMS, PRIME) 
472,360, 'SEED. MAP', 'DEn.HflP', 'PPTflBLE. SEED', 'PRIME' 
m***  STOP 

NOTE: 

The  DEM. MAP  can  be  the  original  DEM,  the  FILLSNGL  processed 
DEM,  or  the  CONCA  VE processed  DEM. 

OK,  SLIST  PPTRBLE.SEED 

1  THIS  IS  THE  COMPLETE  LIST  OF  POUR  POINTS 


SEQUENCE  NO. 

LRBEL 

PRIR 

ELEURTION 

LINE 

SRMP 

1 

0/ 

1 

2070 

2 

0 

2 

0/ 

2 

2151 

4 

0 

3 

0/ 

3 

2146 

7 

0 

4 

0/ 

4 

2126 

7 

0 

163 

55/ 

56 

2608 

448 

49 

164 

0/ 

50 

2485 

453 

0 

165 

0/ 

56 

2381 

469 

0 

166 

0/ 

54 

2194 

463 

0 
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RSNAP  -  Finds  a  cell  with  the  nearest  flow  accumulation  value 
above  a  user  specified  flow  accumulation  threshold  and  within 
a  threshold  distance  from  a  manually  generated  seed  cell.  The 
user  specifies  a  threshold  distance  as  a  value  for  the  number  of 
cells  in  the  line  (row)  and  sample  (column)  direction  that 
represent  the  neighborhood  around  the  seed  cell.  The  program 
then  snaps  a  manually  delineated  seed  cell  to  the  location  of 
the  cell  it  found. 


When  manually  defining  watershed  pour  points  a  problem 
might  arise  in  that  the  newly  created  seeds  are  not  registered  to 
the  modeled  drainage  patterns.  For  example,  a  pour  point 
generated  only  on  the  basis  of  DLG  stream  data  might  fall 
slightly  outside  of  the  drainage  pattern  derived  using  the  HTAS. 
In  such  a  case  the  flow  accumulation  value  of  the  pour  point  is 
low  because  it  is  considered  to  lie  outside  of  the  drainage 
network.  Since  few  cells  are  modeled  as  flowing  into  it,  the 
program  delineates  only  a  very  small  watershed  for  this  pour 
point.  RSNAP  snaps  the  user  defined  pour  points  to  a  cell  on 
the  HTAS  modeled  drainage  pattern  so  that  they  are  properly 
registered  to  the  DEM,  allowing  full  delineation  of  watersheds. 
RSNAP  reads  from  and  writes  to  the  same  file,  therefore  make  a 
copy  of  the  original  file  in  the  HTAS  format.  RSNAP  can  process 
only  up  to  999  points. 


Input  parameters: 

Number  of  lines  or  rows  (NL), 

Number  of  samples  or  columns  (NS), 

The  lowest  flow  accumulation  to  which  the  seed  can  be 
snapped  (FLOW  THRESHOLD) 

The  distance  threshold  from  the  original  seed 
(DISTANCE  THRESHOLD) 

Name  of  the  user  defined  pour  point  or  seed  file  in  single 
quotes:  (SEED  FILE), 

Name  of  the  flow  accumulation  file  in  single  quotes 
(COUNT  FILE), 

Name  of  the  system  in  single  quotes,  here  it  is  'PRIME' 
(SYSTEM) 
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OK,  R  RSNflP 

ENTER  NL,NS,FLOU  THRESHHOLD, DISTANCE  THRESHHOLD,SEED  FILE, COUNT 

FILE,SVSTEri(UNIX,UnS,PRmE) 
472, 36Q ,100,50, 'HflNUflLSEED . HflP ' , 'COUNT . HflP ' , 'PR  I  HE ' 

MOUED  SEED  FOR  85  FROM  30,    251  TO   29, 

CELLS) 

tlOUED  SEED  FOR  65  FROM  66,    207  TO   67, 

CELLS) 

206  TO   69, 


MOUED  SEED  FOR  66  FROM  69, 

CELLS) 

hOUED  SEED  FOR  88  FROM  82, 

CELLS) 

nOUED  SEED  FOR    67  FROM    83, 

CELLS) 


94  TO        83, 
187  TO        83. 


250(1.41421 
207(1.00000 
205(1.00000 
94(1.00000 
188(1.00000 


nOWED  SEED  FOR  109  FROM  358,  146  TO      359,  146(1.00000 

CELLS) 

nOUED  SEED  FOR  106  FROM  362,  102  TO      363,  102(1.00000 

CELLS) 

nOUED  SEED  FOR  111  FROM  421,  144  TO      421,  145(1.00000 

CELLS) 

TOTAL  OF  22  POINTS  NOUED 

114  SEED  CELLS  HERE  FOUND 
****  STOP 
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SEED  -  Finds  starting  or  "seed"  points  for  automatic  watershed 
generation  at  the  outlets  of  potential  stream  tributaries  on  the 
drainage  pattern  modeled  using  the  HTAS.  SEED  accepts  user 
input  of  the  values  (in  number  of  cells)  for  the  flow  accumu- 
lation threshold  above  which  the  pour  points  of  watersheds 
can  be  found.  It  also  uses  flow  accumulation  (COUNT)  values 
and  DELTA  values.  For  each  cell  where  both  the  flow  accu- 
mulation value  and  the  delta  value  are  greater  than  the  flow 
accumulation  threshold,  it  assigns  to  this  cell  a  unique  posi- 
tive value.  This  means  that  for  a  given  threshold  the  program 
places  a  seed  where  the  DELTA  value  or  a  flow  accumulation 
increment  between  two  cells  is  above  that  threshold  and 
where  the  flow  accumulation  (COUNT)  value  is  above  that 
threshold.  The  first  condition  effectively  places  the  "seed" 
where  there  is  a  juncture  or  tributary  relationship  between 
two  potential  streams,  and  the  amount  of  flow  input  (the 
increment)  is  above  the  threshold.  The  second  condition 
eliminates  any  streams  which  have  a  flow  accumulation  value 
below  the  threshold  as  being  too  insignificant  to  have  "seeds" 
created,  thereby  eliminating  minor  tributaries.  The  cells  that 
fit  these  criteria  represent  a  pour  point  or  a  "seed"  of  an 
automatically  delineated  watershed  and  receive  a  positive 
value.  The  watersheds  developed  from  each  of  these  pour 
points  receives  the  value  of  the  "seed"  when  using  the 
WTRSHED  program.  The  background  cells  in  the  SEED  data 
set  receive  a  value  of  -1. 

Input  parameters: 

Number  of  lines  or  rows  (NL), 

Number  of  samples  or  columns  (NS), 

The  flow  accumulation  threshold  for  pour  points 
(THRESHOLD), 

Do  you  want  seeds  at  edges  of  the  data  set:  1  =  no;  0  =  yes 
(EDGES), 

Name  of  the  seed  file  in  single  quotes:  (SEED  FILE), 

Name  of  the  delta  file  in  single  quotes  (DELTA  FILE), 

Name  of  the  flow  accumulation  file  in  single  quotes 
(COUNT  FILE), 

Name  of  the  system  in  single  quotes,  here  it  is  'PRIME' 
(SYSTEM) 

OK,  R  SEED 
ENTER  NL,NSJHRESHHOLD, EDGES, SEED  FILE, DELTA  FILE,   COUNT 
FILE,SVSTEn<UNIX,UMS,PRmE) 

472,360,3000, 1, 'SEED.HflP',  'DELTA. HflP', 'COUNT. HflP',  'PRIHE' 
20  BRSmS  UERE  SEEDED 

!«c)K*>|c    SJQp 
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WTRSHED  -  Finds  watersheds  for  "seed"  cells  that  represent 
user  specified  or  automatically  assigned  pour  locations. 
Delineation  of  watersheds  requires  both  a  flow  direction  data 
set  and  another  "starter"  or  "seed"  data  set  which  represents 
watershed  pour  locations.  The  starter  data  set  consists  of 
background  values  of  -1  in  which  a  "seed"  cell  or  groups  of 
cells  were  inserted  at  the  outflow  points  of  the  desired  water- 
sheds, with  each  start  cell  or  group  of  cells  having  its  own 
unique  positive  value.  The  "seed"  starter  data  set  is  created 
automatically  using  SEED  or  manually  by  generating  the 
points  in  MOSS,  rasterizing  them  to  MAPS  format,  converting 
them  using  MAPSTOIDIMS  and  CONVERT  to  the  direct  access 
file  format  used  by  the  HTAS,  and  snapping  them  to  the 
HTAS  modeled  drainage  pattern  using  RSNAP.  Using  the  flow 
direction  data  set,  "WTRSHED  than  iteratively  reassigns  the 
background  values  of  cells  in  the  starter  data  set  to  the  value 
of  the  "seed"   cell  or  group  of  cells  into  which  the  back- 
ground cells  flow.  The  program  reads  and  writes  to  the  SEED 
file,  so  make  a  copy  of  the  original. 

Input  parameters: 

Number  of  lines  or  rows  (NL), 

Number  of  samples  or  columns  (NS), 

Name  of  the  flow  direction  file  in  single  quotes  (DIR  FILE), 

The  name  of  the  seed  file  in  single  quotes  (MASK  FILE), 

Name  of  the  system  in  single  quotes,  here  it  is  'PRIME' 
(SYSTEM) 

OK,  R  UTRSHED 
ENTER  NL,NS,DIR  FILE, MASK  FILE,SVSTEM(UNlX,UnS,PRIME> 
472, 369 ,  'S6HL .DIR',  'SGHL . SHED ' , ' PR  I  HE ' 


PASS 

1 

PASS 

2 

PASS 

3 

PASS 

4 

PASS 

5 

PASS 

6 

PASS 

7 

PASS 

8 

PASS 

g 

PASS 

10 

PASS 

11 

PASS 

12 

PASS 

13 

PASS 

14 

*))(«*  STOP 
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Resource  &  Environmental  Modeling  Applications 

vuidi  die  HIAS 


There  are  numerous  potential  resource  and  environmental 
applications  of  the  HTA  software  (HTAS).  The  following 
descriptions  of  applications  are  very  general  and  only 
suggest  possible  uses. 

Soli  epQSIon  prediction 

Users  can  model  sheet  and  rill  soil  erosion  with  GIS  using  the 
current  MOSS/MAPS  capabilities  interfaced  with  the  Revised 
Universal  Soil  Loss  Equation  (RUSLE).  However,  the  CONCAVE 
program  in  HTAS  provides  an  added  capability  to  identify 
potential  sediment  sinks  in  the  landscape,  thus  identifying  areas 
of  potential  deposition.  The  software  is  also  being  tested  as  an 
alternative  method  for  determining  slope  lengths.  Currently, 
documentation  is  being  prepared  that  describes  the  GIS  meth- 
odology for  sheet  and  rill  soil  erosion  modeling  with  the 
RUSLE. 

Runoff  modeling  using  tlie  SCS  Curve  Number  method 

The  SCS  Curve  Number  method  generates  storm  hydrographs 
for  watersheds.  It  interfaces  with  GIS  to  help  determine  runoff. 
The  smaller  the  watershed,  the  more  precise  this  information  is 
likely  to  be.  HTAS  permits  subdivision  of  watersheds  into  small 
subwatersheds,  and  thus  provides  a  smaller  analytic  unit  for  this 
type  of  calculation.  Currently  documentation  is  being  prepared 
that  describes  the  GIS  methodology  for  runoff  modeling. 

Sediment  yield  modeling 

Runoff  information  generated  using  the  SCS  Curve  Number 
method  can  be  used  to  provide  the  R  (runofO  factor  for  the 
Modified  Universal  Soil  Loss  (MUSLE)  equation  to  calculate 
sediment  yield  from  each  subwatershed.  Currently  documenta- 
tion is  being  prepared  that  describes  the  GIS  methodology  for 
sediment  yield  modeling  using  MUSLE. 
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Dam  and  reservoir  design 

Using  both  the  MOSS/MAPS  and  HTAS  capabilities  it  is  possible 
to  delineate  watersheds  for  dams,  and  thus  obtain  information 
useful  in  determining  dam  placement.   By  inserting  information 
on  the  elevation  of  potential  dams  into  a  DEM,  and  therefore 
creating  a  depression  in  the  landscape  which  is  then  filled  using 
CONCAVE,  one  can  model  the  area  and  volume  of  reservoirs 
created  by  the  dam  at  any  level  of  water  up  to  the  elevation  of 
the  spillway.  The  resultant  map  may  illustrate  the  areal  extent 
and  volumes  at  various  water  levels.  Currently  documentation  is 
being  prepared  that  describes  this  GIS  methodology. 

Determining  tlie  location 

of  surface  structures  to  controi  erosion 

Flow  information  provided  by  the  HTAS  can  be  used  to  deter- 
mine the  optimum  locations  for  placement  of  erosion  control 
structures,  such  as  spreaders  and  pits.  This  methodology  has 
not  yet  been  developed. 

Drainage  pattern  deiineation/ 
guiiy  and  cliannel  deveiopment 

The  flow  accumulation  map  permits  modeling  of  drainage 
patterns  on  the  landscape,  which  is  particularly  useful  in  the 
absence  of  other  data.  It  will  also  probably  help  in  identifying 
likely  areas  of  gully  and  channel  development,  depending  on 
the  resolution  of  the  digital  elevation  data.  This  methodology 
has  not  yet  been  developed. 


Stream  flow  modeling 

Being  able  to  obtain  values  for  accumulation  of  flow  along  the 
landscape,  when  used  in  conjunction  with  rainfall,  runoff, 
channel  slope  and  other  data,  should  permit  stream  flow  mod- 
eling and  stream  sediment  carrying  capacity  using  combined 
GIS  and  HTAS  capabilities.  This  methodology  has  not  yet  been 
developed. 


Locating  aquifer  recharge  area 

It  is  known  that  aquifer  recharge  areas  are  generally  located 
in  topographic  depressions.  HTAS  provides  that  information. 
"When  used  with  information  on  soil  infiltration  and  presence  or 
absence  of  an  upper  confining  layer,  it  should  be  possible  to 
model  aquifer  recharge  areas  using  GIS.  This  methodology  has 
not  yet  been  developed. 
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Riparian  area  modeling 

Besides  permitting  determination  of  influx  and  routing  of 
sediment  in  streams,  the  HTAS  permits  determination  of  areas 
where  water  might  pond,  loosing  its  energy  and  thus  deposit 
the  sediment.  This  methodology  has  not  yet  been  developed. 


Applications 
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Many  other  applications  are  possible  and  are  being 
explored.  If  you  have  any  suggestions  or  experiences  in  the 
other  applications,  or  are  interested  in  finding  out  more 
about  the  one  described  above,  please  contact: 

Jacek  Blaszynski 
Physical  Scientist 

U.S.  Department  of  the  Interior 
Bureau  of  Land  Management 
Service  Center 
Branch  of  GIS 
Bldg.  50,  SC  344 
Denver  Federal  Center 
Lakewood,  CD  8D225-0047 

303/236-5263 
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Let  us  know. 


This  manual  is  the  result  of  our  effort  to  improve  both  the  content  and 
appearance  of  Service  Center  documentation.  We  have  made  every  effort 
to  ensure  that  this  manual  is  as  functional  as  possible,  but  we  would  like  you 
to  point  out  any  errors  that  you  discover.  This  version  of  the  HTAS  User's  Guide 
contains  many  changes  and  improvements  that  were  inspired  by  comments 
from  the  user  community. 

1 .  Did  you  have  any  difficulty  understanding  or  using  this  manual? 

I i  No  Diificulty  | |  Some  Problems        LJ  Great  Difficulty 

If  you  had  difficulty,  what  specific  areas  caused  you  problems? 


2.  How  can  we  improve  this  manual?  (please  be  specific) 


3.  How  would  you  rate  this  manual? 

n  Excellent      QGood  CD  Fair 

4.  Optional  Information — 

Name: 

Title: 

Address: 

Phone: 

5.  Other  comments — 


•^x 


Please  mail  your  responses  to 

GIS  User  Support  Section  (SC-344C) 

Bureau  of  Land  Management 

Service  Center 

Building  50,  IDenver  Federal  Center 

Post  Office  Box  25047 

Denver,  Colorado  80225-0047 
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